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ABSTRACT 


A  n ew  approach  to  defining  Inlet  distortion  factors  for  axial -flow 
compressors  Is  presented.  A  summary  of  pa;t  analyses  of  Inlet  distortion, 
along  with  a  brief  review  of  flow-distortion  definition  and  compressor 
performance,  provides  a  background  for  the  proposed  theory.  Inlet  flow- 
distortion  In  terms  of  total-pressure  fluctuations  at  the  compressor  face 
Is  converted  to  vortlclty.  The  effects  of  vortlelty  on  the  Internal 
aerodynamics  of  the  compressor  are  Investigated  with  some  approximate 
calculations  Included  to  Indicate  the  validity  of  this  approach.  Refine¬ 
ments  and  alternate  solutions  to  the  theory  which  could  lead  to  definition 
of  a  more  comprehensive  and  reliable  definition  of  stall -Inducing  inlet 
distortion  are  included. 
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I.  INTRODUCTION 


Engine  stall  has  always  been  a  problem  of  varying  magnitude  in  jet 
powered  aircraft  employing  axial -flow  compressors.  At  present  aircraft 
gas  turbine  performance  Is  largely  limited  by  flow  separation  and  surge 
In  the  compressor.  Generally,  highest  compressor  efficiencies  occur 
near  the  point  of  flow  breakdown  (stall),  and  so  the  design  margin 
between  normal  operation  and  stall  Is  kept  to  a  practical  minimum.  A 
compressor's  ability  to  remain  below  this  strll  limit  is  significantly 
affected  by  the  amount  and  degree  of  distortion  in  its  inlet  flow. 

Inlet  flow  distortion  has  been  called  by  many  names,  all  of  which  basi¬ 
cally  describe  variations  in  total -pressure,  total -temperature,  or 
velocity  over  the  inlet  cross-section.  Previously,  inlet  flow  distor¬ 
tion  has  been  at  a  level  which  the  compressor  could  handle  without  a 
serious  degradation  in  performance.  However,  with  the  greater  speeds 
and  performance  demanded  of  today's  high-specific-energy  turbomachines, 
inlet  distortion  at  the  compressor  face  can  alone  be  great  enough  to 
cause  compressor  stall  and  surge.  One  prominent  source  of  distortion 
has  been  traced  to  shock-wave/boundary-layer  interactions  within  the 
Inlet.  Others  are  steam  ingestion,  pressure  transients  due  to  afterburner 
light-off,  ordnance  discharges,  etc.  The  magnitude  of  the  problem  in 
some  recent  aircraft  designs  has  caused  much  research  to  be  conducted  in 
this  area,  specifically  correlating  maps  of  total-pressure  variations  at 
the  compressor  face  to  engine  surge.  The  results  of  this  work  have 
yielded  an  Immense  amount  of  empirical  data.  However,  reduction  of  these 
data  has  yet  to  provide  a  technique  which  adequately  describes  the  distor¬ 
tion  accommodation  of  a  given  engine-inlet  configuration  over  its  entire 
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range  of  operation.  In  contrast  to  empirical  correlations  yielding  a 
statistical  yes-no  answer  to  the  question  of  whether  or  not  the  engine 
will  stall,  methods  are  needed  which  link  stall  criteria  fdr  a  blade 
element  to  distortion  patterns.  This  paper  presents  a  new  approach  to 
analyzing  already  existing  data  and  provides  evidence  of  a  better 
Insight  Into  the  actual  fluid  mechanics  Involved  In  the  problem. 


II.  STATEMENT  OF  PROBLEM 


A.  AXIAL  COMPRESSOR  PERFORMANCE1 

The  performance  of  an  axial-flow  compressor  such  as  found  on  air¬ 
craft  turbojet  engines  Is  best  described  by  Its  steady-state  performance 
map.  Figure  1  Illustrates  a  typical  compressor  map.  "Pressure  ratio" 

Is  exit  total -pres sure  to  Inlet  total -pressure  and  "corrected  airflow" 

Is  engine  mass-flow  rate  corrected  for  total-temperature  and  total- 
pressure. 


O 

h 
< 
q : 


FIG.  I  COMPRESSOR  PERFORMANCE  MAP 
J85-  GE- 13  [adapted  from  Ref.2] 


hhe  Information  and  examples  In  this  section  were  adapted  from  a 
report  by  Brlmelow  [Ref.  1]. 


12 


4 


The  surge  line  defines  the  maximum  pressure  ratios  which  the  compressor 
can  achieve  throughout  Its  range  of  rotor  speeds.  Operation  above  this 
line  Is  accompanied  by  a  drastic  reduction  In  pressure  ratio.  This  1$ 
called  "surge,"  and  It  Is  caused  by  excessive  aerodynamic  blade  stall. 

The  stability  of  a  system  Is  determined  by  the  slope  of  Its  pressure 
ratio/flow  characteristic.  Consider  steady  operation  at  point  of 
Fig.  2. 


CORRECTED  AIRFLOW 


FIG.  2  COMPRESSOR  /  RECEIVER 

MATCHING  CHARACTERISTICS 

Now  If  the  flow  became  unsteady  and  a  perturbation  were  to  cause  a 
reduction  In  airflow  to  point  Sj  {constant  speed),  the  compressor  would 
be  delivering  an  instantaneous  pressure  greater  than  that  in  the 
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receiver.  Effectively  the  compressor  back-pressure  has  been  lowered, 
and  the  flow  would  Increase  until  a  match  Is  again  reached.  Similarly 
for  an  Increase  In  flow  to  point  S2,  the  Increased  back-pressure  would 
cause  a  return  to  stable  operation  at  point  M-|.  However,  If  the  compres¬ 
sor  Is  matched  on  the  positive  slope  of  the  characteristic,  by  the  same 
reasoning,  unsteady  conditions  exist  and  perturbations  are  aggravated 
to  "surge." 

A  stall  margin,  defined  as  the  vertical  distance  between  the  nominal 
(design)  operating  line  and  the  surge  line,  can  be  thought  of  as  a 
measure  of  the  compressor's  ability  to  accommodate  flow  transients  which 
Instantaneously  load  the  machine.  This  margin  Is  dictated  by  overall 
engine-inlet  design  and  Is  affected  by  many  different  phenomena.  Unfor¬ 
tunately  the  point  of  maximum  efficiency  In  a  compressor  Is  just  prior 
to  stall  or  just  below  the  surge  line  on  the  performance  map,  and  there¬ 
fore  a  trade-off  between  efficiency  and  surge  margin  is  required.  The 
expected  flow  transients  which  the  compressor  must  be  able  to  accommodate 
will  affect  the  design  efficiency. 

Flow  transients  Initiated  within  the  engine  system  are  normally 
caused  by  power  level  (throttle)  movement.  Transients  can  also  be 
Initiated  by  flow  distortion  at  the  compressor  face.  Of  the  various 
flow  distortion  parameters,  Brlmelow  [Ref.  l]  states  that  inlet  total- 
pressure  and  total-temperature  have  been  Identified  as  having  the  greatest 
effect  on  engine  stability,  which  Is  to  a  large  extent  the  compressor's 
ability  to  remain  below  the  surge  line.  Total -pressure  variations  can 


i>e  steady-state  spatial  distortion. 


*  i —  - j _ 

Uiwvaf  jf  my  spatial  uuwui  uivut 


In-phase  pulsations  which  are  spatially  uniform  over  the  compressor  face, 
or  a  combination  of  these.  Low  frequency  spatially  uniform  (planar) 
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pressure  fluctuations  (on  the  order  of  20  CPS)  do  not  affect  the  com¬ 
pressor  but  do  have  a  detrimental  effect  on  the  Inlet  and  engine  control 
systems.  High  frequency  planar  fluctuations  (more  than  one  per  revolu¬ 
tion)  affect  compressor  stability  much  like  time-varying  spatial 
distortion. 

Spatial  distortions  can  be  further  segregated  into  radial  and  cir¬ 
cumferential  Inlet  flow  distortion.  This  separation  Into  circumferential 
and  radial  distortion  was  arbitrary,  somewhat  motivated  by  geometric 
simplicity.  A  radial  total -pressure  distortion  such  as  shown  In  Fig. 

3(a)  produces  very  little  variation  In  Inlet  static  pressure.  The  net 
result  Is  a  low  velocity  In  the  low  pressure  region  and  a  high  velocity 
In  the  high  pressure  region. 


FIG.  3  RADIAL  PRESSURE  DISTORTION 

[adopted  from  Ref.  l] 
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In  the  low  velocity  region  the  blades  will  operate  at  higher  Incidences 
and  higher  total -pressure  ratios.  The  net  result  will  be  a  reduction  In 
the  compressor  stall  margin. 

For  the  circumferential  type  distortion  shown  In  Fig.  4(a),  the 
associated  surge  line  reduction  Is  a  function  of  the  basic  stage-match 
of  the  compressor,  the  axial  spacing  between  blade  and  vane  rows,  and 
the  unsteady  response  of  the  rotor  blades.  Since  the  blade  rows  restrict 
circumferential  redistribution  of  flow,  each  circumferential  segment 
performs  essentially  independently. 


FIG.  4  CIRCUMFERENTIAL  DISTORTION 

[adopted  from  Re?.'] 


All  segments  are  required  to  deliver  the  same  exit  static  pressure  so 
that  in  the  area  of  low  inlet  pressure  there  Is  a  high  pressure  ratio 
and  vice-versa.  Correspondingly,  the  high  pressure  ratio  forces  a  low 


axial  velocity  In  this  area  and  a  high  velocity  In  the  segment  of  high 
Inlet  pressure.  These  variations  In  velocity,  produced  less  than  one 
duct-diameter  upstream  of  the  compressor,  result  In  static  pressure* 
fluctuations  at  the  compressor  face.  This  mechanism  of  varying  veloci¬ 
ties  (Incidences)  can  Induce  surge  when  the  average  pressure  ratio  of 
the  distorted  flow  Is  less  than  the  surge  pressure  ratio  for  undistorted 
Inlet  flow.  The  effects  of  circumferential  pressure  pulses  on  a  TF-30 
engine  are  presented  In  Ref.  3.  It  was  found  here  that  the  amplitude  of 
a  pulse  necessary  to  stall  the  engine  varied  Inversely  with  pulse  dura¬ 
tion  and  that  the  engine  was  most  sensitive  to  distorted  sectors  of  180 
to  240  degrees. 

The  distortions  which  are  encountered  In  actual  compressor  appllca-  . 
tlons  such  as  jet  aircraft  Inlets  are  neither  radial  nor  circumferential 
but  a  combination  of  the  two.  Furthermore  the  magnitude,  location, 
shape,  and  duration  of  these  pressure  variations  change  continually. 

It  Is  virtually  impossible  to  evaluate  precisely  the  effect  that  these 
distortions  have  on  engine  stability  because  there  are  so  many  distor¬ 
tion  patterns  significant  to  each  compressor  surge.  Therefore  it  has 
been  necessary  through  experimental  programs  tc  develop  generalized 
correlations  between  pattern  characteristics  and  corresponding  effects 
on  compressor  stability.  Correlating  the  shape,  location,  spatial 
extent,  and  magnitude  of  the  inlet  distortion  requires  analyzing  many 
radial  and  many  circumferential  patterns  and  relating  the  overall  picture 
to  engine  stability.  From  this  is  generated  a  distortion  factor  which  is 
a  measure  of  the  associated  loss  In  stability-margin  for  that  given 
propulsion  system.  Because  compressor  geometry  and  bTadtng  are  not 

considered  In  arriving  at  this  parameter,  the  result  cannot  be  applied 

1 

to  other  systems.  1 
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For  a  more  detailed  and  complete  analysis  of  this  subject,  consult 
Brlmelow's  excellent  work  In  Ref.  1. 

B.  PRESENT  STATE  OF  INLET  DISTORTION  ANALYSES 

Until  a  few  years  ago,  engine-inlet  compatibility  limits  were  success¬ 
fully  determined  by  evaluating  an  engine's  tolerance  to  steady-state 
distortion.  Normally  this  was  done  with  screens  or  similar  obstructions 
placed  In  the  Inlet  and  the  compressor  face  Instrumented  with  low-response 
pressure  probes.  The  engine  was  operated  under  all  conditions  of  Interest, 
and  by  comparing  pressure  data  at  the  compressor  face  with  the  engine's 
performance,  distortion  factors  were  derived.  It  was  known  that  unsteady 
flow  was  being  produced  by  these  test  devices,  but  no  attempt  was  made  to 
evaluate  Its  effect.  However,  when  this  technique  was  applied  to  the 
TF-30  turbofan  engine  during  the  F-'  1A  flight  test  program,  the  data  did 
not  correlate  well  nor  was  the  method  found  valid  for  other  fan  type 
engines  at  that  time.  An  Improved  steady-state  distortion  factor,  K^, 
was  developed,  but  It  too  fell  short  of  providing  the  needed  total  evalu¬ 
ation  of  Inlet  distortion.  It  was  becoming  more  evident  that  there  was 
another  significant  variable  Involved,  namely  turbulence.  Turbulence  can 
be  described  as  flow  with  random  pressure  fluctuations  and  can  be  quanti¬ 
fied  as  a  fluctuation  in  RMS  total -pressure,  (aPrm$)/(pt)*  With 
Increased  aircraft  speeds,  turbulence  has  taken  on  greater  importance 
since  there  Is  a  significant  amount  of  turbulent  flow  associated  with 
the  boundary- layer/shock-wave  interaction  within  supersonic  inlets.  It 
was  felt  that  steady-state  distortion  indexes  were  still  cf  value,  but 
a  correlation  with  turbulence  was  necessary.  One  such  empirical  relation 
showing  a  decrease  In  distortion  tolerance,  Kd2»  with  Increasing  turbulent 
energy  Is  depicted  In  Fig.  5. 
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FIG.  5  LOSS  IN  STEADY- STATE  DISTORTION 
TOLERANCE  WITH  TURBULENCE 
[adapted  trow*  Ref.  4  ] 


Turbulence  has  two  primary  effects  on  axial  compressor  performance. 
First,  turbulent  Inlet  flow  causes  a  reduction  In  compressor  surge  line, 
its  sensitivity  being  a  function  of  turbulent  frequency  and  amplitude. 
Second,  a  reduction  In  flow  capacity  accompanies  Inlet  turbulence.  These 
two  phenomena  were  observed  In  tests  performed  by  Brk:elow  and  Plourde 
[Ref.  5]  on  a  TF-30  turbofan  engine.  This  all  led  to  the  development  of 
new  distortion  Indexes  such  as  the  and  parameters  of  Ref  0. 
However,  lack  of  correlation  caused  Burcham  and  Hughes  to  raise  doubt  as 
to  the  usefulness  of  the  conventional  turbulence  definition,  (aPrm$)/(P-|-)* 
Bellman  and  Hughes  [Ref.  7]  found  no  consistent  relationship  at  all 
between  turbulence  factors  and  surge,  and  in  some  cases  surge  was  experi¬ 
enced  as  turbulence  level  was  decreasing.  Therefore  it  was  concluded  here 
that  peak-to-peak  pressure  fluctuations  at  key  areas  of  the  compressor 
face  are  more  Important  In  surge  prediction  than  are  turbulence  ' wels. 
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Also  extensive  experimental  data  have  been  obtained  from  a  J85-GE-13 

I 

turbojet  placed  In  a  supersonic  wind- tunnel t  running  at  Mach  numbers 
close  to  2.5.  By  varying  the  characteristics  of  the  engine  inlet  duct* 
a  desired  level  of  time*  variant  flow  distortion  was  produced  through 

*  .  t 

Internal  shock  Interactions  [Ref.  2].  Calogeras,  Burstadt  and  Coltrin 
found  for  this  085  configuration  that  time-averaged  distortion  could  be 
well  below  a  stall -Inducing  value  while  instantaneous  peaks  exceeded 

•  i  «  ( 

this  value.  In  addition*  severe  pressure  fluctuations  could  result  in 

i  i 

an  Instantaneous  distortion  of  a  type  quite  different  from  the  time- 
averaged  distortion. 

Of  these  few  representative  programs  discussed,  there  are  two 
striking  similarities  worthy  of  note:  (1)  Each  is  dedicated  to  using 
Inlet  flow  distortion  described  by  raw  total -pressure  data  to  predict 
compressor  surge.  No  attempt,  at  determining  the  cause  or  nature  of  the  ■ 
stall  is  made.  (2)  There  Is  much  uncertainty  and  dlsagreemertt  associated 
with  the  present  methods  of  predicting  distortion  induced  stall.  A  1 

j 

universal  distortion  parameter  has  not  been  agreed  upon. 

I  % 

Ob  •'ously  a  better  understanding  of  the  internal  fluid  dynamics,  1 
while  a  difficult  task  at  best,  would  nonetheless  allow  a  more  compre-  , 
hensive  gei  “ration  of  a  stall-inducing  distortion  parameter.  Total- 
pressure  Is  the  easiest  and  perhaps  t;*e  only  flow  parameter  which  cap  be 
measured  at  the  compressor  face  with  sufficient  accuracy  and  frequency 
response  to  gain  meaningful  data.  The  vast  amount  of  pressure  data 
available  due  to  the  previously  mentioned  programs  provides  an  exten¬ 
sive  empirical  data  base  with  which  to  treat  this  problem.  In 

i 

evaluating  these  data,  inclusion  of  time  and  temperature  dependence 

i  » 

.  I 

could  only  Improve  the  usefulness  of  the  results.  Therefore  It  is 
proposed  that  vorticity  become  the  basis  of  a  "bntversal  distortion 
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factor.1'  This  thesis  will  show  how  existing  total -pressure  maps  along 
with  total -temperature  and  time  data  can  be  transformed  Into  vortlcity 
maps,  Analysis  of  the  resulting  vortlcity  maps  will  provide  evidence 

i 

as  to  the  flexibility  of  this  mejthod. 

Vortlcity  Is  more  closely  related  to  fluid  dynamics  than  Is  total 
pressure  alone.  Compressor  stall  and  engine  surge  are  flow-separation 

i 

phenomena;  therefore  the  key  to  understanding  their  causes  lies  In  the 
Internal  dynamical  situation. 

C.  VORTICITY  APPROACH  ’ 

1.  Compressor  Stall 

.  , 

1  The  most  significant  phenomenon  affecting  the  performance  of  an 
axial-flow  compressor  Is  stall.  If  the  flow  rate  through  a  compressor 

J 

blade  row  Is  decreased  while  the  rotor  speed  Vs  maintained  constant* 

j 

angles  of  attack  on  the  blades  will: increase  and  eventually  flow  sepa- 

,  •  i 

ration-will  occur.  This  condition  is  called  stall  and  is  similar  in 
some  respects  to  the  stall  of  an  Isolated  airfoil;  differences  are  due 
to  the  adjacent  blading  In  a  cascade.  Stalled  operation  of  a  compressor 

1  j  i 

is  accompanied  by  a  drastic  decrease  in  efficiency  and  excessive  oscil- 

i 

lating  blade  loads  capable  of  causing  structural  failure.  Stall  occurs 

i 

In  the  range  of  unstable  compressor  operation  and  hence  can  trigger 
surge.  Usually  a  compressor  stalls  In  an  asynmetric  manner;  that  is, 
well  defined  regions  of  the  annulbs  are  stalled  while  the  remaining  area 
Is  unstalled;  see  Fig.  6.  The  flow  through  these  stalled  patches  is 
severely1 restricted  In  accordance  with  tho  decrease  in  total  flow-rate 
through  the  compressor,  thus  the  flow  Is  diverted  around  this  blockage. 
Effectively  this  Increases  blade  angle-of-attack  on  one  side  of  the 

1  1  i 

restriction  while  decreasing,  It  on  the  other  side.  In  Fig.  7  blade  No.  1 

i 
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i.  6  TYPICAL  STALL  PATTERNS 


Decreasing 
Angle-  of  -  Attack 


Increasing 
Angle-  of  -  Attack 


FIG.  7  STALL  CELL  GROWTH  AND 
ROTATION  IN  A  BLADE  ROW 
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•will  soon  become  unstalled  while  at  blade  No.  5  flow  separation  will 
commence.  This  condition  is  termed  “rotating  stall”  since  it  propagates 
circumferentially  in  a  direction  opposite  to  rotor  rotation.  In  an 
absolute  frame  of  reference  this  speed  of  circumferential  propagation 
Is  30-60  per  cent  of  the  rotor  speed  [Ref.  8],  When  these  conditions 
become  severe  enough  to  cause  net  flow-rate  fluctuations  with  time,  the 
situation  is  called  ''surge.'*  This  describes  operation  in  the  unsteady 
area  of  the  compressor  performance  map  in  Fig.  1.  Valensi  [Ref.  9]  has 
made  a  thorough  investigation  into  these  phenomena  with  some  excellent 
photographic  results.  Of  the  various  analytical  methods  available  for 
predicting  compressor  stall,  the  most  generally  accepted  employs  the 
limiting-biade-loading  criteria  of  Lieblein,  Schwenk,  and  Broderick 
[Ref.  10].  Their  "limiting  diffusion  factor,"  D,  is  a  function  of  the 
velocity  relative  to  the  blade  and  the  circulation  about  the  element. 

If  a  constant  axial  velocity  assumption  is  made,  this  simplifies  to  a 
function  involving  the  flow  turning  angle  and  cascade  solidity.  For 
design  purposes  where  such  stage  parameters  may  not  be  available, 

Lieblein  [Ref.  11]  has  derived  an  equivalent  diffusion  ratio,  D  ,  which 
Is  based  on  inlet  and  outlet  conditions  to  the  compressor  and  the  type 
of  blading  installed.  Unlike  a  single  airfoil  where  angle-of- attack  is 
used  to  predict  stall,  the  critical  variable  for  a  cascade  of  blades  is 
the  flow  turning-angle.  The  exit-flow  angle  as  well  as  the  inlet  angle 
must  now  be  considered. 

2.  Secondary  Flows  in  Axial  Compressors 

Returning  now  to  normal  compressor  operation,  flow  patterns 
based  on  axisymmetric  flow  distribution  and  an  infinite  number  of  blades 
have  been  described  quite  well  in  the  literature.  However,  it  is  necessary 
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to  consider  a  finite  number  of  blades  since  there  are  wall  boundary- 
layer  interactions,  tip  clearance  leakages,  blade  boundary-layers,  and 
centrifugal  effects  which  are  not  represented  adequately  by  assuming 
axi symmetric  flow.  Although  these  asymmetric  phenomena  have  been 
studied  independently  for  reasons  of  simplicity,  it  must  be  remembered 
that  their  effects  are  not  necessarily  additive  (non-linear)  and  may 
have  strong  interactions. 

Figure  8  illustrates  the  various  secondary  flows  which  have 
been  found  to  exist  in  axial -flow  turbomachines  as  a  result  of  the  above 
mentioned  viscous  and  three-dimensional  effects.  Lakshminarayana  and 
Horlock  [Ref.  12]  have  made  athorough  review  of  secondary  flows  in  axial- 
flow  turbomachines  and  have  indicated  where  in  the  literature  more 
specific  information  can  be  found.  Woods  [Ref.  13]  presents  the  current 
state  of  the  art  in  analytical  treatment  of  secondary  flows  and  associ¬ 
ated  losses  in  axial -flow  turbomachines. 

The  induced  velocities  caused  by  secondary  flows  are  considered 
losses  since  only  a  small  part  of  their  kinetic  energy  is  recoverable 
downstream.  Secondary  velocities  may  perturb  the  exit  flow  angle  from 
a  blade  row.  Near  annulus  walls  where  lies  the  origin  of  the  central 
vortex  In  Fig.  8,  the  Induced  exit  flow  angle  can  cause  local  blade 
stall.  Much  work  has  been  performed  in  analyzing  these  phenomena,  and 
some  analytical  methods  for  predicting  three-dimensional  exit  flow 
angles  have  been  proved  valid  through  experimental  results  [Refs.  14 

w.n/4  1 1^1 

WU«  i v j. 

Another  source  of  strong  secondary  flows  in  blade  passages  lies 
In  non-uniform  inlet  flow  such  as  described  previously.  Treating  inlet 
flow  distortion  as  vorticity  provides  a  natural  basis  from  which  to 
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FIG.  8  SECONDARY  FLOW  AND  VORTICES  IN  AN  AXIAL 
FLOW  COMPRESSOR  ROTOR  [adapted  from  Ref.  II  ] 

compute  resultant  secondary  circulation  at  the  exit  plane.  The  effects 
of  this  secondary  flow  on  cascade  performance  can  be  studied  analogously 
to  boundary- layer  induced  circulation.  Comparison  with  an  appropriate 
blade-loading  criterion  may  predict  the  compressor's  unstable  operation. 

4  , 

3.  The  Theory  of  Vorticity  Analysis 

The  thermodynamic  properties  which  are  significant  in  compressor 
performance  are  total -pressure  (velocity),  total -temperature  (enthalpy), 
specific  heat,  and  .molecular  weight.  Variation  of  any  one  or  combination 
of  these  properties  with  time  and/or  position  at  the  compressor  face 
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constitutes  Inlet  flow  distortion  and  should  be  an  Integral  part  of  any 
universal  distortion  parameter.  To  date  distortion  factors  dealing  only 
with  one  Inlet  property  at  a  time  have  been  reported,  l.e.  total- 
temperature  variation  +  Kj  Kj(sta11 )  or  total -pressure  fluctuation  ■* 

Kp  •+  Kp(stall).  This  type  of  piecemeal  approach,  while  providing  satis¬ 
factory  results  for  certain  restricted  applications,  cannot  be  expected 
to  remain  valid  under  all  conditions.  In  reality,  an  aircraft  engine 
will  encounter  all  of  these  fluid  distortion  properties  at  one  time  or 
another.  A  partial  list  of  causes  of  these  Is:  inlet  characteristics, 
environmental  conditions,  rocket  exhaust  from  ordnance,  hot  runway  opera¬ 
tion,  steam  catapult  launches,  and  the  normal  operating  envelope  of  the 
aircraft.  There  Is  definitely  need  for  a  more  realistic  or  universal 
distortion  parameter,  one  which  accounts  for  the  total  inlet  distortion. 

It  is  the  intent  of  this  work  to  show  that  a  distortion  parameter  based 
on  vortlcity  would  be  more  flexible  and  yield  more  comprehensive  surge 
prediction  criteria  than  present  methods. 

Crocco's  Theorem  provides  a  versatile  means  of  converting  exist¬ 
ing  inlet  distortion  maps  to  vortlcity  In  that  pressure,  temperature, 
and  time  are  all  available  as  desired  for  inclusion  in  the  transformation. 
As  indicated  previously,  time  is  an  important  parameter  in  such  an  analy¬ 
sis,  the  aspects  of  steady-state  versus  turbulence  induced  stalls  being 
quite  distinct.  Hawthorne  [Refs.  16  and  17],  Preston  [Ref.  18],  and 
Horlock  [Ref.  15]  all  deal  with  the  passage  of  vortlcity  through  a  cas¬ 
cade.  Her lock  and  Hawthnrnp  calculate  Induced  exit  flow  angles  which 
Involves  a  solution  to  Poisson's  Equation.  Calculating  appropriate  exit 
circulation  and  comparing  it  with  blade-ltft  circulation  is  a  far  simpler 
process  that  yields  a  qualitative  feel  for  the  relative  size  of  this 
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effect.  Comparison  of  this  secondary  circulation  with  that  caused  by 
three-dimensional  effects  should  also  be  made  In  order  to  look  for 
possible  Interactions. 

The  effect  of  Inlet  distortion  on  blade  angle-of-attack  must 
also  be  considered  since  this  is  the  other  “half"  of  the  flow  turning 
angle.  This  effect  Is  much  less  subtle  and  can  be  dealt  with  In  a  more 
straight-forward  manner.  Areas  of  lew  pressure  represent  lower  axial 
velocity,  hence  higher  angles-of-attack.  The  unsteady  response  of  the 
blades  must  be  considered,  particularly  for  circumferential  Inlet 
distortion  where  blade  passing  frequency  Is  a  factor. 


HI.  CONVERSION  OF  PRESSURE  TO  VORTICITV  HAPS 
A.  CROCCO'S  THEOREM 

In  general  form,  Crocco's  Theorem  can  be  written  as 
T  grad  S  +  u  x  curl  u  *  grad  H  +  ££  (1) 

p  t 

which  holds  In  unsteady  flow  [Ref.  19].  From  basic  Thermodynamics 


Combine  Eqns.  (1)  and  (2),  assuming  H  *  Constant 

-jr-  grad  P-j-  +  u  x  curl  u  *  ^  (3) 

T  3t 

The  following  development  will  apply  to  the  system  defined  In  Fig.  9. 


r 


C  -  Chord 


FIG.  9  NOMENCLATURE 
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The  following  reference  quantities  taken  at  the  compressor  face  will 
be  used  to  non-dimensional ize  the  equations.  Refer  to  Fig.  9. 


r  *  p'Fr 

Vj  =  Average  Stagnation  Pressure 

u  *  u'a 

a  =  Average  Speed  of  Sound 

r  «  r'| 

D  =  Compressor  Diameter 

z  *  z‘  . 

M  ,  i  2W 

“it 

»  =  Fluid  Vorticity 

Table  I  lists  some  of  the  characteristic  times  available.  Blade  passing 
frequency  will  be  used  in  this  development,  since  most  available  data 
Include  compressor  speed;  then 

*  -  *■  •  2ttR 
v  w  V  N 

tip" 

The  numerical  value  of  the  reduced  frequency  for  stall  is  based  on  the 
work  of  Carta  [Ref.  20].  For  values  less  than  approximately  0.1,  blade 
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TABLE  I 

CHARACTERISTIC  Tl  MES 

BLADE  PASSING  FREQUENCY 

vt.pN  . 
2.w  R 

(1500)  (36) 

(6)  (2)  '  4500  Hz 

ONCE-PER-REV 

V—1D 

1500 

-J.LP-  - 
2trR 

mz)  "25Hz 

WAVE  PROPAGATION 

RADIAL  — 

2R 

1200 

■  4  »,  300  Hz 

LONGITUDINAL 

a  +  u 

— 8  60,180 

«  120  Hz 

STALL  TIME 

REDUCED 

»  JZl£  4  0.1 

FREQUENCY 

u 

.  L  (0.1)  (1200) 
(5)  {0.2/ 

8  200  Hz 
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stall  Is  essentially  equivalent  to  steady-state  stall.  For  a  reduced 
frequency  larger  than  0.1,  stall  Is  suppressed  due  to  dynamic  effects. 
Other  work  In  this  area  has  been  performed  by  Carta  [Rpf.  22].  The 

i 

paper  by  Ericsson  and  Reding  [Ref.  21]  Is  a  survey  of  dynamic  stall. 

Rewriting  the  component  equations  In  non-dimensional  form  and 
canceling  terms  yields 


1  3p' 
yP'  3r' 


VnpNa 


(7) 


1  ap' 
yPV  ae 


+  -=•  (uz'u>r'  -  u'«  ') 
a 


vTIPNa 

2ira2 


(8) 


1  aP' 
yP'  az‘ 


r 


U6V> 


VTIPWa  /3“z'\ 
2ua^  \  3t '  / 


(9) 


Equations  (7),  (8),  and  (9)  are  exact.  Orders  of  magnitude  of  various 
terms  will  be  examined  to  ascertain  whether  or  not  complexity  can  be 
reduced  in  a  logical  and  consistent  manner. 

B.  ORDER  OF  MAGNITUDE  ANALYSIS 

Figure  A-2  represents  ’n  Inlet-flow  distortion  map  at  the  face  of 
a  J85-GE-13  turbojet.  This  distortion  map  triggered  surge.  The  map 
of  Fig.  A-3  is  a  similar  situation  in  a  TF-30  turbofan  engine.  These 

2 

Reference  to  figures  and  equations  labeled  with  a  letter/number 
combination  nave  more  detailed  explanations  which  are  contained  in  the 
Appendix. 
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maps  will  be  considered  representative  of  stall-inducing  flow  distortion 
for  purposes  of  making  an  order-of-magnltude  analysis  on  Eqns.  (7),  (8), 
and  (9). 


Establish  an  order  of  magnitude,  e  ,  such  that 


From  Figs.  A-2  and  A-3 


-0.2  <  e  <  0.2 


Between  12  and  5  o'clock.  Fig.  A-2,  3P'/3r'  Is  a  maximum.  Refer  to 
Appendix  A  for  all  numerical  values. 

3P‘  afPv/F,)  1.17  -  .97 

- -  =  - ! — L  *  -  *  i.53 

3r'  3(2r/D  (7.1  -  6.05)(2/16.1) 


From  Fig.  A-3  at  the  12:30  o'clock  position 


3P'  1.04  -  .96 

-  « - r  .685 

3r'  (2.15  -  1.8) (2/6) 

3P' 

-0.7  <  - <1.5 

~  3r'  “ 

Between  5  and  7  o'clock.  Fig.  A-2,  a  maximum  value  of  aP '/ae  exists  such 

that  1.20  -  .95 

—  s  -  «  .416  per  radian 

ae  .6  rad 


Similarly  from  Fig.  A-3  at  1  o'clock 


aP  *  1.04  -  .92 

—  - - =  .343  per  radian 

36  .35  rad 


.4 


\_ 


36 
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An  average  value  of  uz*  for  the  J85  can  be  obtained  from  'data  in 
Appendix  C. 

u,  460 

u  '  =  JL  “ -  *  .415 

2  a  1110 

A  maximum  value  would  be  i 

1/2 

*  750  ft/sec 


Data  with  which  to  estimate  the  magnitudes  of  ur*  and  ue'  or  their, 
partial  derivatives  are  not  available.  However,  an  upper  bound  of  order 
c  can  be  definitely  established  by  considering  the  consequences  of  ur’ 
or  ue'  =  0(e°).  A  radial  velocity  component  o.f  order  unity  is  not 
possible  within  a  duct.  Any  ur*  components  present  at  the  compressor 
face  are  caused  solely  by  turbulent  fluctuations  and  are  probably  of 
order  e^;  however  an  upper  bound  of  e  will  be  assumed.  A  velocity  ■ 
component  such  as  u0',  called  pre-rotation,  is  often  generated  intention 

i 

ally  as  a  means  of  establishing  a  desired  angle-of-attack  on  the  first 
blade  row.  The  magnitude  of  this  velocity  should  be  smaller  than  uz' 
and  certainly  not  larger.  A  velocity  in  the  circumferential  direction 
of  order  e®,  about  equal  to  the  axial  velocity,  Is  unrealistic;  and  If 
tt  were  present;  the  engine  would  just  not  run.  For  the  caseiunder 
consideration  there  is  no  pre-rotation  so  that  any  u0'  which  does  exist . 

o 

is  caused  also  by  turbulent  fluctuations  and  Is  probably  of  order  zu,. 
Again  an  upper  bound  of  e  will  be  assumed.  1 


uz  3 

(PT  *  PS) 
P/2 

750 

u2*  * - 

=  .675 

1110 

uz'  -  .7 
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Therefore 


Ua'  «  U  '=  <  U  ' 

e  r  z 


Partial  derivatives  of  u0*  and  ur'  with  respect  to  time  can  also  be 


given  only  a  bound  with'respept  to  uz'. 


3u0'  ,  3ur*  3uz' 


Therp  are  two  phenomena  causing  a  time  dependency.  One  is  due  to  the 

i  .  • 

potion  of  the  rotor.  For  patterns  of  large  extent  on  the  Pj  map,  the 

i 

appropriate  frequency  is  once-per-rev.  For  patterns  of  limited  extent, 

i  t  , 

blade  passing  frequency  is  appropriate.  Unsteady  aspects  of  the  flow 
cause  a  time  dependency  also. 

,  l 

Continuing,  the  magnitude  of  wr‘  will  bp  assessed. 
wJ)  /  3u_  3y/\  D 

V  =  —  *  (  +  /&)•-= 

1  2a  i  \  38  Az  /  2a 


From  Fig.  A-2  between  8  and  9  o'clock,  and  data  from  Appendix  C 

(?  -  p  y/2 

P*  =  .97  u  *  '-I - V  =  400  ft/sec 

2  p/2 


P‘  =  1.17 


uz  =  750  ft/sec 


A8  =  30°  *  .52  radians 


(350  ft/sec) (1.34  ft) 

(.52  rad) (.605  ft)2(lil0  ft/sec) 


s  .814  per  radian 


“r  «  .8 


I 
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Similarly,  approximating  yields 
UgO  /  3Up/  3uz  \  D 
we  ~  'll  lr  /  2? 


For  the  axial  component  of  vorticity 

1  &izD  /  3u  3ur  \  D 

«z'  “  —  =  (  —  +  —  )  ” 

2a  \  3r  30  /  2a 

J 

Since  both  partial  derivatives  are  really  unkown  in  size,  this  term 
will  be  evaluated  later. 

The  magnitude  of  3uz'/3t'  in  Eqn.  (9)  r.an  be  estimated  by  again 
referring  to  the  point  "fi"  in  Figs.  A-3  through  A-8. 


u2  *  231  ft/sec 

P*  *  .96  u2  =  374  ft/sec 

At  «  (6  figs.)  x  (.2x10" 3  sec/fig.) 

3u2*  it  /Au2\  3.14  /  3^  ~  231  \ 

7t~  a  RPS  NUt/  (1190)  (166)  \  1.2xl0”3  / 

3uz'  ,  .06 
TT 
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Entropy  can  be  expressed  as 


»»here  P  and  T  are  static  values  at  the  compressor  face.  Assuming  P  =  F, 
combination  of  this  equation  with  Eqn.  (2)  yields 


1  + 


4 


1  + 


Y-l 
(e)  Y 


1  + 


2 

e 


The  deviation  of  local  speed  of  sound  from  average  speed  of  sound  is  less 
than  the  deviation  of  local  stagnation  pressure  from  average  stagnation 
pressure.  Therefore  a/a  will  be  set  equal  to  1. 

Table  II  is  a  summary  of  these  established  relative  mrgnitudes  in 
terms  of  e.  The  velocities,  u0'  and  ur' ,  are  placed  at  thjir  absolute 


upper  bounds  and  wz 

'  is  set  equal  to 

e  for  later  verification. 

TABLE  H 

RELATIVE 

ORDERS  OF 

MAGNITUDE 

€* 

€ 

€* 

€*' 

P' 

d? 

TY 

jlp!  jip; 
ae  TP 

a -a 

a 

u  '  u  ' 
ur  e 

“z' 

GDr*  CUe‘  CDj 

at’  av 

duz' 

at' 
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Now  Eqn.  (7)  can  be  written 


1  3P‘  ,  N 

T7  +  (ueV  -  W)  * 


yP*  3r* 
e°  e° 


0  0 
e  e 


VTIP  N  3ur‘ 
2w  a  3t* 


Dropping  terms  of  order  e 


J_  3P|_ 

Yp*  3r* 


“e  = 


«2V  e  0 


— f-) 

J,u  •  \  ar‘ / 


Similarly  for  Eqn.  (8) 


1  3  P 1 

yP'r*  ae 


+  (uz'u>  '  -  U-*  w7 ' )  = 


c°e° 


VTIP  w  3V 
2 if  a  at* 

0  2 

e  e 


1  /3P'\ 

V  =  — : — f— )  (11) 

yP,uz'r'\9&  / 

The  size  of  w  '  was  assumed  to  be  of  order  e°  although  at  this  point  it 
is,  1r  fact,  unknown.  In  arriving  at  Eqns.  (10)  and  (11),  ojz *  is 
multiplied  by  small  terms  (uf‘  and  ug');  therefore  little  if  any  error 
should  be  introduced  here.  Subsequent  analysis  establishes  the  size  of 
&>2'.  From  Eqn.  (9) 


1  aP'  . 

”  77T  +  (ur'V 

It  -A.  * 

i  I  I  I 


Un  V 


VTIP  ^  auz* 
2n  a  at1 


I  I 


e°  « 
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Equation  (12)  contains  most  of  the  terms  whose  magnitudes  were  generally 
uncertain.  All  terms  are  of  the  order  e  so  no  consistent  simplification 
can  be  made.  This  equation  will  not  be  used  in  the  subsequent  trans¬ 
formations  from  Pj  to  vortlclty. 

C.  APPLICATION  OF  RESULTS 

Equations  (10)  and  (11)  are  the  required  relations  for  transforming 
total -pressure  maps  to  circumferential  and  radial  vortlclty  maps 
respectively.  The  procedure,  as  applied  to  the  J85  map,  is  outlined 
below. 

The  expressions  3P ' /9r;  in  Eqn.  (10)  and  3P'/3e  in  Eqn.  (11)  are 
obtained  in  a  similar  fashion  from  Fig.  A-2.  For  one  particular  radial 
line,  values  of  P'  are  plotted  against  r',  the  slope  at  any  point  being 
the  local  value  of  ap'/ar'.  This  is  done  for  sufficient  radials  to 
adequately  cover  the  entire  map.  Local  values  of  P*  can  be  taken  direct¬ 
ly  from  the  Py  map,  and  from  these  uz'  can  be  calculated  (assuming 
constant  static  pressure).  Similarly  in  Eqn.  (11),  for  one  discrete 
radius,  values  of  P'  at  a  position,  0,  are  taken  at  Intervals  through 
360  degrees.  These  values  are  plotted,  the  local  slopes  being  numerical¬ 
ly  equal  to  3P'/a0.  Other  values  necessary  to  compute  «r*  can  be  taken 
directly  from  the  map  or  readily  computed.  See  Appendix  B  for  the 
calculations  outlined  above. 

The  results  of  these  calculations,  depicted  in  Figs.  10  and  11  as 

and  ujg *  contour  maps,  present  a  two-dimensional  picture  of  conditions 
at  the  compressor  face.  Since  v*(vxu)  *  0  ,  it  follows  from  Gauss' 
Divergence  Theorem  that 

->• 

//  («  •  r)  dA  *  in  (v  •  oj)  dV  *  0 
A  y 
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Application  of  this  theorem  to  the  maps  of  &>r'  and  contours  yields 
a  quantitative  estimate  of  the  axial  vortlcity,  uz',  which  Is  reflected 
In  Table  II.  These  calculations  are  contained  also  in  Appendix  B,  and 
the  resulting  map  is  given  as  Fig.  12. 

Radial  and  circumferential  vortlcity  maps  for  the  TF-30  total- 
pressure  data  in  Fig.  A-10  have  also  been  generated.  They  are  Figs. 

13  and  14. 
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FJG.  M  CIRCUMFERENTIAL  VORTICITY  CONTOURS,  J85 
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IV.  ANALYSIS  OF  VORTICITY  MAPS 


A.  RADIAL  VORTICITY  INLET  EFFECTS 

In  vector  notation  radial  vortlcity  is  aligned  parallel  to  local 
blade  spans  of  the  compressor,  +wr  pointing  outward  radially  according 
to  the  right-hand  rule.  The  vorticity  in  Fig.  10  is  concentrated  in 
two  areas,  -uir  between  5  and  7  o'clock  and  +wr  at  8  o'clock.  In 
analyzing  the  effects  of  these  local  concentrations  it  must  be  noted 
that  four  distinct  situations  can  be  represented  by  a  vorticity  map  such 
as  this.  These  are  listed  in  Table  III.  Vortex-sheets  are  aligned  par¬ 
allel  with  the  flow  and  vortex  filaments  are  convected  with  the  flow. 


TABLE  III 

Possible  Combinations  for  Radial 

Inlet  Vorticity 

Case 

Vortex  Blades 

Flow 

I 

Sheet 

Stator 

Steady 

II 

Sheet 

Rotor 

Unsteady 

III 

Filament 

Stator 

Unsteady 

IV 

Filament 

Rotor 

Unsteady 

Case  I 

A  continuous  sheet  of  vortex  filaments  passing  through  a  stationary 
blade  row  (stator)  Is  a  steady-stat-.  problem.  According  to  the  Dict- 
Savart  Law,  velocities  are  induced  in  the  fluid  which  cause  changes  in 
local  blade  incidences.  Fig.  15.  In  a  cascade  r  is  the  circulation 
about  a  blade. 
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FIG.  15  RADIAL  VORTEX  "SHEET  PASSING 
THROUGH  A  STATOR 


Therefore  blades  in  region  "A"  where  ar/r  Is  positive  are  more  heavily 
loaded. 

The  first  stator  row  is  behind  a  rotor  so  Figs.  10  and  13  cannot  be 
applied  here  without  knowing  the  influence  of  the  rotor  on  u>r'.  This 
case  could  be  applied  directly  to  inlet  guide-vanes  however. 

Case  II 

A  steady  vortex-sheet  passing  through  a  rotating  blade  row  is  an 
unsteady  problem  due  to  the  relative  motion  between  particle  paths  and 
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For  opr  sign  convention;  causes  blade  unloading  while  +wr  Increases 

i 

the  local  Incidence.  Comparison  of  circulation  causing  blade  lift  to 
that  associated  with  Inlet  distortion  '(rdlst)  ^rom  ^9«  W  should  give 
h  feel  for  the  relative  strength  of  this  effect.  For  brevity,  all  mean 

i 

quantities  will  be  used  and  the  dynamic  response  of  the  blading  will  be 
neglected.  From  data  in  Appendices  C  and  D,  for  the  J85 

i 

Wui  =  850  ft/sec 

'  1 

ku2  =  321  ft/sec' 

Blade  Gap  (Sm)  =  ft 

i 

i  Blade  Chord  (C)  *  — —  ft 

1  12 

'  riift  ■'  S(Wul  -  Wu2)  '  (16) 

I  I 

rii.ft  =  3°  ft2/sec  (Relative) 

» 

rd1st  =  “rcs  08) 

1) 

From  Fig.  A-2  at  8  o'clock 


2854.9 _ 

P'AUP'  -8.17 


1  ,  3P 


-(-) 
r'  '  36  • 
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The  ratio  of  circulations  becomes 
r11ft  30 

- -  *  -  a  10 

rdist  2.97 

which  is  about  10  per  cent.  It  must  be  remembered  that  this  Is  ;curring 
at  blade  passing  frequency,  and  any  Induced  stall  would  be  dynamic  stall. 

Dixie  Passing  Frequency  -  ^ TSTT  7,3  TST- 

Case  III 

Unsteady  vortex  filament  motion  through  .  stator  is  a  more  complex 
problem  since  time  dependency  of  the  fluid  motion  as  well  as  the  unsteady 
response  of  the  blading  must  be  considered.  This  situation  can  be 
pictured  as  the  pass? je  of  a  vortex  filament  Jr  a  local  concentration  of 
vorticity  passing  through  a  cascade.  In  a  sim'lar  problem,  Crimi  and 
Reeves  [Ref.  25]  calculated  sizable  lift  fluctuations  as  a  helicopter 
rotor  blade  passed  over  a  spanwise  oriented  vortex.  Ward  [Ref.  26] 
observed  dynamic  stall  of  a  rotor  blade  which  he  asserted  was  caused  by 
a  radial  vortex  passing  beneath  the  blade.  Figure  1/  depicts  this 
situation  in  a  cascade. 


FIG.  17 


RADIAL  VORTEX  FILAMENTS 
FASSING  THROUGH  A  STATOR 


Case  IV 


It  follows  that  a  vortex  filament  approaching  a  rotating  blade  row 
Is  also  an  unsteady  problem.  Figure  18  shows  basically  how  Wul  is 
affected  by  wr. 


RADIAL  VORTEX  FILAMENT 
APPROACHING  A  ROTOR 


B.  CIRCUMFERENTIAL  VORTICITY  INLET  EFFECTS 

Circumferential  vorticity  taken  vectorially  is  positive  in  the 
clorkwise  direction.  Due  to  the  signs  in  Eqn.  10,  the  positive  gradients 
in  Figs.  11  and  14  represent  -juig'J-  The  circumferential  vorticity  map 
for  the  J85  in  Fig.  11  shows  a  strong  concentration  of  vorticity  at  the 


outer 


- - -14. 

CA  LC11U  1  I 


f Fiv^minh  about 

Wilt  VMJJM  -  - 


'inn  dearees. 


Figure  14  shows  a 


similar  characteristic  for  the  TF-30  although  it  is  weaker  and  less 
extensive.  As  in  the  case  of  „v\  four  basic  situations  can  exist  for 
possible  analysis.  These  are  listed  in  Table  IV. 


TABLE  IV 

Possible  Combinations  for  Circumferential 

Inlet  Vorticity 

Case 

Vortex  Blades 

Flow 

1 

Sheet 

Stator 

Steady 

II 

Sheet 

Rotor 

Steady 

III 

Filament 

Stator 

Unsteady 

IV 

Filament 

Rotor 

Unsteady 

Case  I 

Circumferential  vorticity  is  oriented  perpendicular  to  the  local 
blade  spans  of  tne  compressor.  A  sheet  of  u0  approaching  a  stationary 
blade  row  is  pictured  in  Fig.  19. 


FIG.  19  CIRCUMFERENTIAL  VORTEX-  SHEET 
APPROACHING  A  STATOR 
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Induced  velocities  according  to  the  Law  of  Biot-Savart  cause  Va^  and 
hence  a-j  to  change  ahead  of  the  cascade.  This  results  in  variable  lift 
along  the  blade  spans  which  produces  shed  vorticity  downstream.  Shed 
vorticity  will  be  dealt  with  In  more  detail  later  when  considering 
secondary  flows. 
v^se  II 

A  circumferential  vortex  sheet  approaching  a  rotating  blade  row 
presents  two  problems  for  analysis.  Relative  motion  between  the  particle 
paths  and  the  rotor  gives  rise  to  a  streamwise  component  of  vorticity  as 
well  as  a  component  perpendicular  to  the  blade  spans.  Fig.  20. 


VORTEX- SHEET  APPROACHING  A  ROTOR 

The  concentration  of  w  '  at  the  outer  radius  in  Fig.  11  is  essentially 
constant  for  about  300  degrees  of  angular  displacement.  The  life  time 
of  one  particular  characteristic  such  as  this  can  be  as  long  as  the  time 


for  one  revolution.  The  series  of  TF-30  Pj  maps  in  Appendix  A  shows  a 
life-time  of  close  to  1/4  revolution.  The  significance  of  this  is  that 
one  blade  can  experience  the  same  velocity  distortion  for  up  to  2*50 
degrees  of  revolution  which  allows  a  steady-state  analysis.  Induced 
velocities  caused  by  the  normal  component  of  this  sheet  are  those  shown 
in  Fig.  19.  Velocities  induced  by  the  streamwise  component  are  pictured 
in  Fig.  21. 


FIG.  21  STREAMWISE  COMPONENT  OF  A 

CIRCUMFERENTIAL  VORTEX- SHEET 
AT  A  ROTOR 


The  streamwise  component  causes  the  typical  variations  in  lift  along  the 
blade  span,  but  it  also  passes  through  the  blade  row  and  affects  the 
exit  flow  angle.  This  will  be  discussed  later.  The  normal  component 
also  causes  a  spanwise  variation  in  lift  as  shown  in  Fig.  19. 

For  of  small  circumferential  extent,  the  situation  is  unsteady 
due  to  relative  motion  between  particle  paths  and  the  rotor. 
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Case  III 


A  circumferential  vortex  filament  approaching  a  stator  is  an  un¬ 
steady  problem.  The  effect  of  the  associated  entry  velocity  perturbations 
shewn  in  Fig.  22,  a  function  of  the  dynamic  response  of  the  blading,  would 
probably  be  negligible.  However,  as  this  filament  passes  through  the 
blade  passage,  a  streamwise  component  develops  which  can  interact  with 
existing  secondary  flows.  This  phenomenon  will  be  discussed  in  the  next 
section. 


FIG.  22  CIRCUMFERENTIAL  VORTEX  FILAMENT 
APPROACHING  A  STATOR 


Case  IV 

A  circumferential  vortex  filament  approaching  a  rotor  is  different 
than  the  previous  case  only  because  a  streamwise  component  of  the  filament 
is  present  ahead  of  the  blades.  Again  unsteady  response  of  the  blades  is 
a  factor  to  consider. 
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C.  PASSAGE  AND  EXIT  EFFECTS 


Several  authors  have  treated  the  problem  of  turning  a  non-uniform 
flow  in  a  cascade  of  airfoils.  Squire  and  Winter's  method  for  calcu-. 
lating  secondary  flows  downstream  of  a  cascade  [Ref.  27]  has  been  used 
extensively  in  the  literature.  Preston  [Ref.  18],  using  a  more  simpli¬ 
fied  approach,  also  treates  the  generation  of  streamwise  vorticity  which 
gives  rise  to  secondary  flow  after  a  cascade.  However,  Preston  assumes 
that  in  addition  to  the  components  of  vorticity  present  in  the  original 
vortex  filament,  a  change  in  circulation  about  the  blade  will  cause 
circulation  to  be  shed  into  the  stream.  Hawthorne  [Ref.  16  and  17] 
attempts  an  analytical  proof  of  the  existence  of  these  different  com¬ 
ponents  of  vorticity  in  nonuniform  flow,  and  he  obtains  reasonable 
agreement  with  experimental  results. 

Consider  an  Isolated  vortex  filament  in  a  sheet  normal  to  the  flow 
direction  as  it  passes  through  a  cascade  of  stationary  airfoils;  see 
Fig,  23.  Due  to  the  difference  in  flow  velocity  over  the  suction  and 
pressure  surfaces  of  the  blades,  this  filament  will  be  turned  and 
stretched  as  It  passes  through  the  cascade.  Downstream  of  the  cascade, 
segment  ab'  gives  rise  to  a  "distributed  streamwise  vorticity"  through¬ 
out  the  passage  exit  plane.  This  description  applies  not  only  to  a 
filament  as  shown  here  but  also  to  a  vortex-sheet.  Squire  and  Winter's 
well  known  expression  for  this  distributed  vorticity  is 


’dist 


=  -2 z&L 
dz 


(19) 


where  £  is  streamwise  vorticity,  e  is  the  deflection  angle,  and  Z  is 
the  spanwise  position. 


FIG.  23  VORTEX  FILAMENT  PASSAGE 
THROUGH  A  CASCADE 


Segment  b-b',  part  of  the  original  filament  a-b  in  Fig.  23,  lies  along 
the  stagnation  streamline  and  in  considering  the  spanwise  direction,  this 
becomes  a  "vortex-sheet."  This  is  called  "trailing  filament"  vorticity. 

A  third  component  of  exit  vorticity,  "shed  circulation"  caused  by  vari¬ 
able  spanwise  lift  due  to  the  asymmetric  inlet  flow,  also  lies  along  the 
stagnation  streamline  and  gives  strength  to  the  trailing  vortex-sheet. 

This  component  Is  labeled  c-d  in  Fig.  23. 

Preston's  simple  analytical  expressions  for  these  secondary  components 
of  circulation,  as  applied  to  Fig.  23,  are  as  follows. 

Distributed  Secondary  Circulation: 

dr-,  =  £?  S  cos  (£.)  dz 

<  -  z 


And  replacing  with  Squire  and  Winter's  expression  for  distributed 
secondary  vorticity,  Eqn.  (19) ,  yields 

drl  =  ~2e(~dzOS  C0S  ^  dz 


(20) 
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r  iiifn-ii  u 


Trailing  Filament  Circulation: 


dr2  sin  (-|_)  dz  (21 

Trailing  Shed  Circulation: 

dr3  *  ( 2  S  sin  (_|_)  dz  (22) 

Trailing  shed  circulation  is  equal  in  magnitude  and  direction  to  dr2. 

The  net  circulation  associated  with  the  trailing  vorticity  in  one  blade 
passage  and  from  one  airfoil  in  the  form  of  a  "vortex  sheet"  is  the  sum 
of  these  three  components. 

dr  -  dr^  +  dr2  +  dr3 

dr  =  Gii1)  C  4  s  s1n  (f)  -  2  e  s  cos  (f)  3  dz  (23 

From  this  exit  secondary  circulation  the  change  in  mean  exit  flow  angle 
can  be  calculated.  This  total  circulation  is  not  uniformly  distributed 
over  the  exit  area  [Ref.  18]. 

Investigations  into  changes  in  deflection  due  to  cascade  secondary 
flows  have  generally  taken  two  lines  of  thought.  Values  of  induced 
deflection  at  the  lifting  line  were  derived  from  the  Induced  velocities 
calculated  solely  from  the  strength  of  the  vortex  sheets  by  Carter 
[Ref.  28]  and  Hausmann  [Ref.  29]  to  mention  a  few.  These  methods  neglect 
the  effects  due  to  distributed  passage  vorticity  and  hence  are  unlikely 
to  provide  satisfactory  agreement  with  experiment  [Ref.  12], 

Hawthorne  [Ref.  30],  Smith  [Ref.  31],  and  others  have  considered 
secondary  motions  to  be  induced  by  the  distributed  vortices  within  the 
passages.  It  Is  assumed  that  the  flow  downstream  behaves  as  if  each 
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blade  passage  were  separated  from  its  neighbor  by  the  blade  wakes,  i.e. 
vortex  sheets.  The  secondary  velocities  in  the  blade  passages,  v 
circumferentially  and  w  spanwise,  are  defined  by  a  stream  function,  if>, 
such  that 


3^ 

w  =  "  3y 


3Z 


The  stream  function  satisfied  continuity  for  this  steady  secondary  motion 
superposed  on  the  main  flow. 


3w  +  iv  =  o 
ay  az 


Substituting  the  stream  function  into  the  vorticity  equation  (u>  =  curl  u) 
and  setting  it  equal  to  Squire  and  Winter's  expression  for  the  distributed 
vorticity,  Eqn.  19,  yields  Poisson's  Equation. 


a2* 

- —  T  — 5* 

ay*  az* 


a2ip 


■  -5*  =  -2e^ 


(24) 


Solutions  to  Eqn.  (24)  are  tedious  but  have  been  accomplished  by  a  series 
method  [Refs.  27  and  30],  by  relaxation  techniques  [Ref.  17],  and  by 
electrical  analogs  [Ref.  32].  The  latter  reference  considers  variable 
spanwise  deflection  which  causes  Eqn.  (24)  to  become 


V2*  -  *  >,$  (25) 


An  exit  flow  solution  for  a  parabolic  inlet  profile  similar  to  the  case  % 

at  hand  is  given  by  Honda  [Ref.  33].  % 

* 

i 

i 

£ 

l 

f. 

\ 


Once  secondary  velocities,  v  and  w,  have  been  determined,  local 
changes  in  the  exit  flow  angle,  Aa2 ,  can  be  computed.  V2  is  the  main¬ 
stream  exit  velocity  for  uniform  inlet  flow  with  exit  angle  a2. 


The  exit  flow  angle  for  distorted  inlet  flow  then  becomes 

a2  «  a2  +  a«2 

dist 

Comparison  of  a2  with  appropriate  blade  loading  criteria  [Ref.  10 
dist 

and  11]  could  explain  the  compressor's  unstable  operation. 

Another  method  of  dealing  with  secondary  velocities,  v  and  w,  ob¬ 
tained  from  solving  Poisson's  Equation  was  used  by  Hawthorne  [Ref.  17]. 
Again  the  assumption  of  blade  passages  extended  downstream  by  vortex- 
sheets  trailing  the  blades  is  made.  There  are  no  normal  velocities  at 
the  vortex-sheets.  The  strength  of  these  sheets  can  be  determined  by 
calculating  aw,  the  change  in  spanwise  secondary  velocity  across  them. 

As  mentioned  before,  these  sheets  of  fluid  leaving  the  trailing  edges 
are  vortex-sheets  composed  of  trailing  shed  and  trailing  filament 
vorticity;  see  Fig.  24.  If  the  trailing  filament  vorticity  is  estimated 
and  subtracted  from  the  total  strength  of  the  sheet,  the  trailing  shed  , 
circulation  remains.  From  this  can  be  calculated  the  variation  in  lift 
on  the  blade  since  shed  circulation  and  blade  lift  are  related  by 

Al  =  Pur3 

Hawthorne's  expression  for  this  trailing  filament  circulation  is 

dr2  ■  -.»1J  *i  .  -2£cos(.s)l  (26) 

C  \  dz  /j_  2  COS  aj  *•  J 

\ 

\ 

\  * 
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FIG.  24 


u ROTOR 

SECONDARY  FLOW  DOWNSTREAM  OF  A 
BLADE  ROW  [  adapted  from  Ref.  17] 


which  Is  the  same  as  Preston's  formula  given  as  Eqn.  (21)  except  for 
different  inlet/exit  flow  angles.  Hawthorne's  predicted  shed  circulation 
agreed  quite  well  with  an  experimentally  determined  spanwise  lift 
variation. 

D.  CALCULATION  OF  SECONDARY  CIRCULATION 

A  calculation  of  these  components  of  secondary  circulation  as  caused 
by  the  circumferential  vorticity  in  Fig.  11  was  performed  in  order  to 
quantitatively  evaluate  the  effect  of  such  distortion  on  a  compressor. 

The  analysis  was  made  assuming  steady-state  conditions  which  means  Case 
II  of  Table  IV,  a  circumferential  vortex-sheet  passing  through  a  rotor. 
The  following  terms  require  definition. 

V  -  Absolute  Velocity 
W  -  Relative  Velocity 
U  -  Rotor  Peripheral  Velocity 


w  -  Absolute  Fluid  Vortlclty 
C  -  Relative  Fluid  Vortlclty 
n  -  Rotor  Angular  Speed 
a  -  Absolute  Flow  Angle 
3  -  Relative  Flow  Angle  , 

-*•  -  Velocity  Vector 
■*+  -  Vortlclty  Vector 
Subscripts;  see  Fig.  25 

r,$,z  -  Absolute  Coordinate  Direction 
s,n  -  Natural  or  Streamline  Coordinates 
s  -  Streamwise 
n  -  Normal  to  Streamlines 
0  -  Upstream  of  Compressor 

1  -  Compressor  face  (entry  to  blade  row) 

2  -  Blade-row  exit  plane 


\ 
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■  '  I 

Horlock's  method  [Ref.  15]  of  treating  the  passage  of  a  vortex 
filament  through  a  rotating  blade-row  results  in  the  following  vector 
transformation.  Figure  26(a)  is  the  velocity  diagram  ahead  of  the  rotor 
which  is  convectlng  the  vortex  filament,  w0  (one  of  many  In  a  sheet). 

The  absolute  velocity  Is  axial  so  that  W0-|  =  Figure  26(c)  is 

the  Inlet  absolute  vorticity  combined  with  rotor  rotation  to  give  the 

i  > 

inlet  relative  vorticity.  Equations  (27)  and  (28)  are  the  corresponding 

—  .  , 

I 

analytical  expressions  for  these  components  of  relative  vorticity  in 

i 

terms  of  mfl. 

I 

i 

1 

CS1 1  =  +2  ft  cos  3j  -  ui0  sin  3^  (27) 

tnl  *  -2  ft  sin  31  -  o>0  cos  3]  (28) 

The  approach  of  a  vortex  filament  to  a  rotating  cascade  and  the  resultant 

i  • 

secondary  relative  vorticity  at  the  passage  exit  are  pictured  in  Fig.  27. 
This  figure  is  similar  to  Fig.  23  except  it  applies  to  a  rotor  while  Fig. 
23  is  for  a  cascade. 

Combining  inlet  streamwise  vorticity  which  passes  through  the  cascade 

i 

unaffected  with  the  exit  vorticity  of  Eqn.  (19)  yields  the  total  dis¬ 
tributed  vorticity  in  the  blade  passage  exit  plane. 

i 

5s2  =  Csl  "  2  e  Si  (29) 

Substituting  Eqns.  (27)  and  (28) 

i 

Cs2  =  (20.  cos  3)  -  u>e  sin  3])  -  2e(-u0  cos  3]  -  2n  sin  3])  (30) 

Since  this  exit  vorticity  is  in  the  direction  of  the  flow,  multiplying 
it  by  the  blade  passage  exit  area  normal  to  the  direction  of  flow  yields 

I 

exit  distributed  circulation. 

i  > 
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FIG.  27  PASSAGE  OF  A  VORTEX  FILAMENT 
THROUGH  A  ROTATING  CASCADE 


df~  **  ?s2Scos62  (31) 

When  transformed  from  absolute  to  relative  coordinates,  Hawthorne's 
expression  for  trailing  filament  circulation  given  previously  as  Eqn. 
(26)  becomes 


sin  232  -  sin  2b^ 
cos  0i 


2e  cos  02J  (32) 


Assuming  two-dimensional  flow,  Hawthorne's  equation  for  shed  circulation 
is 


dr 


5  cos  »i  <^i  ( tan  -  tan  S}) 


\  ■**  ~*  / 


A  simplified  application  of  these  equations  to  the  circumferential 
vorticity  and  total -pressure  maps  for  the  J85  (Figs.  11  and  A-2)  was 


made  by  assuming  an  average  radial  distribution  of  3P'/3r'  and  P'  to 
exist  for  a  circumferential  segment  of  the  face.  Equation  D-2  was 
applied  to  radials  at  the  three  circumferential  positions  shown  in 
Fig.  28.  (Appendix  D) 

»e  -  (  -2854.9  \  3P1  sec->  (D-2) 

\  PV9.17P'  -  8.17/  ar 

Now  using  Eqns.  (27),  (28),  and  (29),  the  various  components  of  relative 
vorttclty  were  calculated  for  a  sufficient  number  of  radial  position* 
see  Appendix  E.  Figure  29  shows  the  relative  sizes  of  these  components 
for  circumferential  position  I. 


Secondary  circulations  as  a  function  of  span  at  the  blade  row  exit  plane 
were  calculated  using  Eqns.  (31),  (32),  and  (33);  see  Appendix  E.  The 
area  under  curve  "A"  in  Fig.  30  represents  distributed  secondary  circu¬ 
lation.  The  area  under  curve  "B"  is  trailing  filament  circulation,  and 
the  area  under  "C"  is  shed  circulation.  These  latter  two  are  the  bir- 
culations  due  to  the  vortex-sheet.  The  magnitudes  of  these  components 
are  found  by  integration  over  dr'. 

.9 

R  /  / flLiA  dr*  =  32.1  ft2/sec 
.3  \  dr  / 

where  R  is  the  radius  used  for  non-dimensional i zing. 

Concerning  the  results  in  Table  V,  distributed  secondary  circulation, 
r-j,  represents  a  loss  in  kinetic  energy.  Blade-lift  circulation  is 
approximately  equal  to  30  ft2/sec;  therefore  the  ratio  of  these  two  is 
about  1  to  1.  In  order  to  evaluate  the  actual  effect  r-j  has  on  the 
compressor,  a  solution  to  Poisson's  Equation,  as  described  previously, 
would  be  necessary.  From  this  solution  the  exact  strength  of  the 
trailing  vortex-sheet  could  be  determined.  Subtracting  trailing  fila¬ 
ment  vorticity  (or  from  the  calculated  strength  of  this  vortex-sheet 
would  yield  the  exact  shed  circulation.  Equation  (33)  from  which 
in  Table  V  has  been  calculated  was  found  by  Hawthorne  to  yield  values 
higher  than  those  found  experimentally.  However,  this  analysis  is 
essentially  an  order-of-magnitude  check,  and  Eqn.  (33)  is  adequate. 

Since  shed  circulation  is  related  to  the  spanwise  variation  in  blade 
lift  by 

At  =  p  r3 

AL  „  rshed 
L  r11ft 


then 


(34) 


™  *■  -n ’*  r** **•»  i'> yrof  .r,’tvT!>,~fyy'T~<> 


The  circulation  about  one  blade  in  the  rotor  under  consideration  for 
uniform  inlet  flow  is  approximately 

r11ft  ■  S<“U1  -  wu2>  *  30  ft2/sec 

The  average  value  of  positive  shed  circulation  over  the  three  circumfer¬ 
ential  data  points  considered,  Radials  I,  II,  and  III,  is 


5.0  +  3.9  +  15.1 
3 


8  ft^/sec 


So  the  vocation  in  lift  from  Eqn.  (34)  is 


hi  8 

T  ’  30 


25  % 


which  represents  unloading  at  the  center  span  of  the  blade.  For  the 
negative  shed  circulation  at  the  blade  tip 

-6.0  -  8.5  -  10.9  AL  -8.8 

*  *  »  on  v 


If  the  compressor  blading  is  operating  near  stall  it  certainly  seems 
plausible  that  a  change  of  25  to  30  per  cent  in  lift  could  trigger  stall. 
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V.  SUMMARY  AND  REVIEW 


Figures  31  and  32  are  presented  as  a  review  and  outline  of  the 
preceding  theory  for  analyzing  the  effect  of  inlet  flow  distortion  on 
an  axial-flow  turbomachine.  Data  from  total -pressure  maps,  3Pj/ar  and 
3pj/ae,  consitute  the  initial  known  quantities.  All  pertinent  assump¬ 
tions  and  simplifying  restrictions  are  included  in  the  chart  as  they 
apply.  The  resulting  determination  of  wz  through  a  consistent  order- 
of-magnitude  analysis  implies  complete  definition  of  the  velocity  field 
ahead  of  the  compressor  with  the  exception  of  the  ambiguity  in  sign  of 
ioz.  This  is  significant  in  that  a  three-dimensional  flow  field  has  been 
determined  from  a  two-dimensional  pressure  map  at  the  expense  of  intro¬ 
ducing  an  error  of  order  e2.  Continuing  in  Fig.  32,  the  analysis  next 
relates  portions  of  several  existing  theories  on  secondary  flows  in 
cascades  and  turbomachines  to  the  maps  of  radial  and  circumferential 
vorticity.  This  yields  an  estimate  of  the  strength  of  the  secondary 
flows  caused  by  inlet  flow  distortion.  The  large  value  for  AL/L  which 
is  calculated  suggests  that  the  linear  cascade  theories  used  are  not 
adequate  for  the  accuracy  required  in  deriving  a  stall  criterion.  Some 
necessary  considerations  for  further  refinement  of  the  theory  which  were 
not  considered  in  this  analysis  are  (1)  blade  twist,  (2)  three-dimensional 
effects  and  interaction  between  normal  secondary  flows  with  distortion 
Induced  flows,  and  (3)  changes  in  vorticity  with  passage  through  the 

blade  i*OW. 


FIG.  3!  FLOW  DIAGRAM  FOR  ANALYSIS;  STEADY  FLOW 


Absolute  Vorticity,  Map 


Relative  Vorticity  Components: 

Horlock  [Ref.  15]  51 


?s2  =  fUnl* 

Squire  &  Winter  [Ref.  27] 


Total  Exit  Streamwise  Vorticity,  cs g 


Circulation  Components,  r 
Hawthorne  [Ref.  16] 
Preston  [Ref.  18] 


Poisson’s  Eqn.  -*•  Induced  Velocity 
v  =  Normal  to  Span 
w  *  Spanwise 


Lift  -  f (rshed) 
AL/L  ■  Arshed/r 


Vortex  Sheet  -*  aL 
Hawthorne  [Ref.  17] 


Change  Exit  Flow 
Angle,  a&£ 


Limiting  Blade  Loading 
Criteria,  Lieblein[Refs.  11  &  12] 


UNIVERSAL  DISTORTION  FACTOR 


Fig.  32.  Composite  Theory  for  Analyzing  Yorticity  Maps 
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VI.  CONCLUSIONS 


This  analysis  Is  the  first  attempt  at  relating  total -pressure 
distortion  maps  to  the  Internal  aerodynamics  of  a  turbomachine.  Sepa¬ 
rating  the  Inlet  distortion  Into  components  of  radial  and  circumferential 
vortlclty  Introduced  flexibility  In  treating  the  problem  as  well  as 
maintained  a  comprehensible  physical  picture  of  the  aerodynamical 
situation.  For  the  data  analyzed  herein,  radial  vorticity  was  found 
to  be  less  significant  to  degradation  of  compressor  performance  than  was 
circumferential  vortlclty  since  the  radial  vorticity  was  coupled  to  the 
dynamic  response  of  the  blading.  This  would  require  an  unsteady 
analysis.  The  relative  significance  of  wr  and  in  this  analysis  does 
not  necessarily  apply  in  general.  The  effect  of  axial  vorticity,  <$z. 
must  be  included  in  further  work. 

The  size  of  the  calculated  shed  circulation  from  a  blade  indicates 
distortion  Induced  secondary  flows  are  of  sufficient  size  to  seriously 
affect  the  compressor's  operation.  With  further  refinement  to  the 
present  theory,  secondary  flows  could  be  evaluated  to  a  degree  of 
accuracy  capable  of  predicting  stall  through  the  mechanism  of  shed 
circulation  or  changes  in  exit  flow  angles.  A  solution  to  Poisson's 
Equation  would  provide  the  needed  accuracy  with  which  to  use  appropriate 
blade-loading  criteria  In  defining  inlet-distortion  stall  criteria. 
Although  stall  normally  originates  somewhere  downstream  of  the  first  few 
stages,  the  distortion  must  oe  treated  as  it  passes  through  each  indi¬ 
vidual  stage.  The  present  theory  deals  only  with  the  first  blade  row; 
however,  it  Is  one  row  closer  to  the  origin  of  the  stall  and  has  elimi¬ 
nated  a  certain  degree  of  uncertainty.  This  theory  may  be  extended  to 
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succeeding  stages  in  the  compressor  by  the  technique  illustrated  in 

i 

Fig.  26.  ' 

Finally,  in  working  toward  a  universal  distortion  factor,  the 
effects  of  other  types  of  inlet  distortion,  i.e.  total -temperature  and 
molecular  weight,  must  be  included.  These  extensions  to  the  present 
theory  should  provide  reliable  and  comprehensive  stall -inducing 
distortion  criteria. 


APPENDIX  A 


i  ORDER  OF  MAGNITUDE  ANALYSIS 

: 

i  i  1 

Part_J_  Preparation  of  Total -Pressure  Data 

The  following  data  and  calculations  apply  to  the  experiments  per¬ 
formed  on  a  085-GE-13  turbojet  as  described  in  Ref.  2,  a  =  0°  stall 
point.  Some  data  were  also  obtained  from  Ref.  24. 

Subscripts:  J  0  -  Free  stream 

1  -  Compressor  face 

2  -  Rotor  exit,  first  stage 
T  -  Stagnation  quantity 

i  Tj]  =,1549  psf 

Rj  =  138C  psf  (Static) 

1  1  .  IT1  =  532° R 

'  M  =  2.5 

:  co 

2 

Compressor  Inlet  flow  area,  A-j  =  185  in 
Figure  A-l  is  a  pressure  contour  map  of  Pt^/Tsto  reproduced  from 
Ref.  2.  For  our  purposes,  this  must  be  converted  to  P'  which  is  defined 
as  PT1/Ptl  . 

\ 

Pressure  Recovery  of  Inlet  =  .77 


.77 


FIG.  A-l  PEAK  INSTANTANEOUS  PRESSURE 
CONTOURS,  J  85  ( from  Rof.  2  ) 


FW"3V’TT'i,1F; 


Application  of  Eqn.  (A-l)  to  Fig.  A-l  yields  the  desired  pressure 
contours.  Fig.  A-2.  Figures  A-3  through  A-l 1  are  instantaneous  pressure 
contours  (P ‘ )  taken  from  experiments  conducted  on  a  TF-30  engine  [Ref. 
23].  The  radial  lines  at  the  centers  represent  compressor  rotation. 

RPM  *  10,000 

time  between  Figs.  =  / \/  _L  £§¥.  \ 

\104  Rev /\  32  seg  / 

=  .1875  x  10“3 


=  .2  ms. 


Elapsed  time  (A-3  ->  A-12)  -  1.5  ms. 

Figure  A-ll  is  the  frame  prior  to  engine  surge. 

F$  =  1440  psf 

Estimated  Inlet  Values 

FT  =  .4 

Y 

PT  =  Ps  (  1  +  41  M2)  Y_1 

Ff  =  1610  psf 
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Flow  Into  Paper 


FIG.  A- 2  TRANSFORMED  PRESSURE 

CONTOURS 
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FIG.  A-6  FRAME  FOUR  OF  INSTANTANEOUS 

PRESSURE  CONTOURS  SEQUENCE,  TF-30 
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F!G.  A- 7 


FRAME  FIVE  OF  INSTANTANEOUS 
PRESSURE  CONTOURS  SEQUENCE,  TF  -  30 


i _ . 
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Part  II  Pressure  Fluctuations  in  the  Z  Direction 

In  the  case  of  steady  flow,  APy  is  present  only  when  dissipation 
occurs.  Some  causes  of  such  dissipation  in  an  inlet  are  shock  waves 
flow  separation,  obstructions,  and  boundary  layers.  These  examples 
are  illustrated  in  Fig.  A-12.  _ 


Boundary 
Layer  _ 


/ 


Flow 

Separation 


(b) 


Wall  Effects 


oundary  Layer 


Strut 


(c)  Obstructions 


FIG.  A- 12  SOURCES  OF  DISSIPATION 

WITHIN  •  INLETS 

♦ 
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The  stagnation  pressure  gradient  is  present  only  in  the  vicinity  of  dis¬ 
sipation.  Therefore  in  steady  flow,  3 P .' / 3 z ‘  is  zero  at  the  compressor 
face  if  the  dissipation  is  located  upstream  as  shown  in  Fig.  A-13. 


t 


FIG.  A-13  PRESSURE  HISTORY  FOR 
STEADY  INLET  FLOW 

For  unsteady  flow,  pressure  gradients  are’  again  caused  by  dissipation. 
If  there  were  a  standing  shock  wave  in  the  inlet,  unsteady  flow  could 
cause  this  to  oscillate  in  the  axial  direction.  This  is  called  "shock 
wave  twitter"  and, can  be  self  excited  in  oscillation.  This  is  depicted 
in  Fig.  A-14. 

Xo 


Range  of 


FIG.  A-14  SHOCK- WAVE  TWITTER 
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The  changing  geometry  seen  by  this  oscillating  shock-wave  causes  pres¬ 
sure  pulses  to  be  generated  as  shown  in  Fig.  A-15.  A/A*  is  the  area 
ratio  as  a  function  of  x.  The  second  quadrant  is  local  Mach  number  as 
a  function  of  area  ratio,  A/A*.  The  third  quadrant  is  the  stagnation 
pressure  ratio  for  a  shock  at  M.  Finally,  the  fourth  quadrant  yields 
the  derivative  of  PT  with  respect  to  x. 


FIG.  A-15  TOTAL  PRESSURE  FLUCTUATIONS 
FROM  OSCILLATING  SHOCK -WAVES 

"  i 

The  resulting  aPj/ax,  a  function  of  diffuser  design,  can  be  evaluated 
at  t Jie  compressor  face  by  assuming  an  acoustic  analysis  with  no  wave 
reflections  at  the  compressor.  Figure  A-16  is  a  plot  of  the  pressure 
pulses  (waves)  in  the  x-t  plane.  The  slope  of  the  characteristics  is 
equal  to  the  reciprocal  wave  speed;  curvature  is  due  to  the  variable 
area  at  the  center-body. 
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(time) 


FIG,  A-16 


PRESSURE  WAVES  IN  THE  x-f  PLANE 
AS  A  RESULT  OF  SHOCK- WAVE  TWITTER 


Consider  the  three  points  labeled  1,  2,  and  3  in  Fig,  A-16. 


z2  "  2]’  (z2  -  )(2/D)  1  (z2  -  z.) 

V  ■  V  ^3  ■  a  (t3  -  t}) 


1 

=  (a  +  u)  =  a'  +  u' 
a 


therefore 


and 


!!1 

3Z' 

3PT 


_J _ /  3P'\ 

t'  +  U1)  l  3t '  J 


a*. 


3z  (a  +  u)  \  3t 


(A-2) 


The  usefulness  of  Eqn.  (A-2)  can  be  demonstrated  by  application  to  a 
simple  inlet/engine  configuration  with  an  oscillating  shock  wave  present; 
see  Fig.  A-17.  For  brevity  the  following  assumptions  are  made 
(1.)  a  =  Constant 

(2.)  Wave  impedance  of  compressor  matches  wave  impedance 
of  duct  v 

(3.)  U  =  constant 

(4.)  "No  wave  dispersion 


k - L  - 


FIG.  A-17  NOMENCLATURE 
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The  shock  wave,  oscillating  within  ±N,  produces  pressure  waves  which 
are  represented  as  a  sine  wave  function. 


x  -  x0  a  N  sin  (a>t) 

s  111  S  COS  (ust) 

within 

gPT  ax  3Pt 
it”  "  9t  ix” 

*  as  N  COS  (list) 

xo 

Using  a  wave  transit  time  equal  to  L/(a+u)  with  a  «  wave  speed  relative 
to  the  gas,  the  unsteady  pressure  at  the  compressor  face  becomes 


ifi 

3t 


=  U  N  COS  U)(t  - 


a?T 

W 


*0 


Combining  this  with  Eqn.  (A-2)  yields 


Non-dimensional i zed  the  preceding  equation  becomes 


*"•  >  / 
where  n  Is  rotor  angular  speed. 

This  Is  one  method  of  relating  pressure  fluctuations  at  the  compressor 
face  (aP'/az’)  with  the  pressure  waves  generated  by  Inlet/shock-wave 
interaction. 


aP' 
a?  t 


(w/n)(flN/a)  cos  [  fl  ^  +  a*  +  u* 


Lfl 

I")] 


( 
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APPENDIX  B 

CONVERSION  OF  Pr  MAP  TO  VORTICITY  -  J85 

Pressure  data  taken  directly  from  Fig.  B-l  are  listed  In  Table  B-I. 
From  these  data  the  curves  in  Fig.  B-2  are  generated.  Slopes  taken 
from  these  curves,  which  are  tabulated  in  Table  B-II,  provide  a  means 
for  constructing  the  circumferential  vortlclty  map  In  Fig.  11.  A 
similar  procedure  is  performed  in  generating  the  radial  vortlclty  map 
in  Fig.  10.  Data  are  contained  in  Tables  B-I II  and  B-1V  and  Fig.  B-3. 
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Table  B-I  Pressure  Data  Taken  In  Radial  Direction  from  J85 


Total -Pressure  Map 


0‘ 

r' 

> 

P' 

30 

r’ 

0 

P! 

60 

r‘ 

0 

P' 

90( 

r' 

> 

P‘ 

120° 

r  P‘ 

150° 

r‘  P' 

.4 

.91 

.47 

.91 

.46 

.91 

.44 

.91 

.71 

.91 

.47 

.84 

CO 

V£7 

• 

.97 

.7 

.97 

.75 

.97 

.62 

.91 

.75 

.97 

.55 

.84 

.72 

1.04 

.74 

1.04 

.79 

1.04 

.72 

.91 

.79 

1.04 

.68 

.91 

.78 

1.10 

.79 

1.10 

.82 

1.10 

.75 

.97 

.82 

1.10 

.75 

.97 

.81 

1.17 

.81 

1.17 

.84 

i  .17 

.79 

1.04 

.84 

1.17 

.78 

1.04 

.92 

1.17 

.92 

1.17 

.92 

1.17 

.82 

1.10 

.92 

1.17 

.80 

1.10 

.95 

1.04 

.95 

1.04 

.95 

1.04 

.85 

1.17 

.95 

1.04 

.82 

1  17 

.92 

1.17 

.92 

1.17 

.95 

1.04 

.95 

1.04 

98 


Table  B-II  Pressure  Gradients  Obtained  from  Slopes  in  Fig!  B-2 


r* 

.4 

.5 

.6 

.7 

.76 

.78 

.80 

.82 

.84 

.87 

.9 


.94 


.96 


0° 

ap’/ar' 

30° 

r'  ap'/ar' 

60° 

r'  3P'/3r' 

90° 

i 

r'  3P'/3r' 

120° 

r'  3P'/3r’ 

150° 

r*  3P73r 

0 

.4 

.13 

.4  .07 

.32 

.4  .03 

.3! 

-.22 

.09 

.5 

.16 

.5  .12 

.38 

.17 

.5  .06 

\ 

.3 

-.22 

.28 

.6 

.26 

.6  .135 

•  .44 

.26 

.6  .1 

.5 

i 

-.22 

.9 

.7 

.5 

.7  .25 

.62 

m 

.68  .33 

.51 

1.8 

.76 

1.4 

.74  .56 

,6 

-.26 

.72  1.08 

.54 

.3 

2.6 

.79 

5 

.8  2.3 

.62 

-.27 

.8  3 

.6 

.53 

1.6 

.8 

6 

.82  4.5 

.66 

D 

.82  2.5 

.7 

.73 

.83 

.82 

1.1 

.86  .55 

.7 

.29 

.84  1.08 

.76 

1.5 

.37 

.86 

0 

-.89  0 

.72 

1:2 

.06  .4 

1 

.79 

4.0 

0 

.9 

-.72 

.92  -.6 

1 

.75 

2.0 

!  .89  0 

.8 

7 

-.6 

.94 

-2.7 

.94  -2.7 

.8 

2.0 

.9  -.26 

8.2 

1.3 

-2.7 

.96 

-3 

.96  -3.0 

.83 

6.0 

.92  -.6 

.86 

0 

-3.0 

..85 

.9 

r 

.94  -2.7 

.9 

r.7  2 

.94  -2.7 


.96  -3.0 


Table 'B- 1 1  (Continued) 


180°  ! 

r'  aP'/ar1 

!  210°  , 

r'  aP'/^r* 

240° 

r'  aPVar' 

270° 

r'  aP */3r’ 

300° 

r'  aP '/ar ' 

330° 

r'  aP'/ar' 

.34 

-.25 

.3  ’  .07 

.3 

.13 

.32 

.505 

.4  .35 

.4  .15 

.4 

-«35 

.3  .07 

i .  • 

.3 

.13 

1 

.4 

.38 

.5  .26 

.6  .33 

148 

-.12 

1.0 

| 

1 

.5 

.28 

.6  .35 

.7  .88 

.51 

0  1 

> 

:  .9 

.13 

.6 

.35 

.7  .55 

.76  1.9 

.6 

.37 

.92 

0 

.7 

.52 

.74  .8 

.78  2.1 

.7 

.7  . 

1 

.94  - 

.15 

1 

.76 

1.6 

.78  2.3 

.8  1.08 

.78 

1.7 

i  i 

1  • 

H 

.78 

.56 

.8  1.3 

.82  .6 

1 

|. 

l 

Hgfig 

.15 

.8 

.44 

.82  .65 

.84  .23 

.86 

1.7 

.85 

0 

.85  0 

.86  0 

.88, 

EX 

.89 

-.37 

.9  -.7 

.9  -.72 

.9 

-2.5 

i 

.92 

-.6 

.92  -1.3 

.94  -2.7 

.96 

-1.2 

1 

.94 

-2.7 

.94  -2.7 

.96  -3. 

i 

.96 

-3.0 

.96  -3. 

\ 


Table  B-III  Pressure  Data  Taken  in  Circumferential  Direction  from  J85 

Total -Pressure  Map 


•  r'  =  .46  r*  =  .54  r‘  =  .61 


.69  I  r*  =  .77 
P' 


4.19  .91  1.81  .91  1.69  .91  0  .9 

4.48  .97  2.35  .84  2.41  .84  1.25  .91 

5.23  .97  3.14  .84  3.19  .84  2.79  .91 

.174  .91  4.10  .91  3.98  .91  3.84  .91 

4.41  .97  4.36  .97 

5.48  .97  5.72  .97  4.59  1.04 

5.13  1.04 


4.48  1.04 


0  P' 
(rad) 


— 

— 

.174 

1.17 

3.70 

v91 

.785 

1.10 

3.92 

.91 

1.39 

1.04 

4.22 

.97 

2.0 

1.10 

4.39 

1.04 

2.26 

1.17 

4.52 

1.10 

2.61 

1.17 

4.71 

1.17 

2.87 

1.10 

4.88 

1.17 

3.14 

1.04 

i 

5.23 

I 

1-17 

3.40 

.97  6.11 

1.17 

r'  =  .93 

0  P' 
(rad) 


3.03  1.17 
3.19  1.10 
3.37  1.04 
3.58  .97 

4.10  .97 

4.27  1.04 
4.45  1.10 


.It 


Table  B-IV  Pressure  Gradients  Obtained  from  Slopes  in  Fig.  B-3 
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Application  of  Gauss'  Divergence  Theorem  Is  performed  In  the  follow¬ 
ing  manner.  A  grid  such  as  Is  shown  In  Fig.  B-4  Is  placed  over  a 
vortlclty  map. 


FIG.  B-4  METHOD  OF  DIVERGENCE 

APPLICATION 


The  net  u>r'  passing  through  region  abed  Is  equal  to 


[(-.I3)(2fe)  +  (2^!2)(i&]  -  [(J£0)(!  cd)  +  (2iiia)(3  cd)]  .  «_• 


Net  flux  of  t»e'  Is  calculated  In  a  similar  manner  using  the  same  grid 
over  an  «6*  map.  The  sum  of  flux  of  <*>r'  and  «0'  is  then  equal  to 
w2'(ab‘Cd)  so  that 

uz'(ab‘cd) 

(It  '  a  - - 

z  Area l ahrd 1 

»  *  -  -  -•  * 


This  procedure  Is  performed  on  each  of  the  36  segments  so  that  a  complete 
wz'  map  can  be  drawn.  These  data  are  contained  In  Table  B-V.  If  the  sum 
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of  net  flux  of  «r’  and  w0‘  is  negative,  then  the  net  flux  of  «2'  roust 
be  positive.  There  Is  an  ambiguity  of  algebraic  sign  which  cannot  be 
resolved  without  additional  data  or  Information.  One  evaluates  . 
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Table  B-V  Net  Fluxes  of  Vortlcl 


+.022  0  -.39  36  0  -.003  +.05 


APPENDIX  C 


FIG.  C-l  BLADE  PASSAGE  DIMENSIONS 

Calculated  Parameters 

5a  .  (,♦*!*,  A 

jfgj  =  1.125  -  (1  +  .2M2)3'6 
^  =  .414 
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TTZ  *»  R  a" 


a  *  49.02  y5U 


a  *  1110  ft./sec 


IL  =  ( .414) (1110  ft/sec) 


U,  =  460  ft/sec 


Pj  * 


p  = _ im.jw.ft2. _ 

(1715.6  )  (514°R) 

d  =  .001565  slugs/ft3 


As  a  check  on  these  figures,  another  method  may  be  used. 


jjcor  <  „  (33.5  lbs/sec)  \  2TT6  ' 


r532/518 


W  -  24.15  lbs/sec 


W  *  pA }iiz 


_ \~.zi  r.."/ _ _ v.:../. a 

(.001565  slu^s  )  (155  in2)  (32.2 


»  462  ft/sec 
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Estimation  of  Sj 


FIG.  C-2  VELOCITIES  AHEAD 

OF  ROTOR 


Estimation  of  $2 

Individual  stage  characteristics  and  loss  data  for  the  J85  are  not 

available.  The  following  data  are  from  Ref.  34: 

Compressor  of  Free-Vortex  Design 

Llebleln  Diffusion  Parameter*  D  <  .5 

Solidity  (a)  =  1.14  -i 

y  First  Stage  Rotor 
Axial  Velocity  Ratio  ®  .96  J 

D  -  (1  -  Hi  )  +  (Ref.  10) 

■1  2YW} 

W2  =  Va2/cos  e  *96  Val/cos  Bp 

*  442/cos  Bp  »  546 

Wu2  *  Va2  tan  e2  *  .96  Val  tan  Bp 

*  442  tan  Bp 

Wi  *  Ya/cos  P-j  *  802  ft/sec 

5  B  a  _  442/cos  Bp  )  +  (650  -  442  tan  Bz) 

802  (2) (1 . 14) (802) 

Bp  B  '36°  minimum 

Wu2  «  320  ft/sec 

Since  the  loading  limit  Is  D  <  .5,  36°  can  be  considered  a  minimum  for 
B2»  *nd  the  maximum  turning  angle  becomes 
AS  «  p  -  Bp  =  55w  -  36° 

-  19° 

e  «  .332  radians 

t 

i 
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Substituting  into  E;?n.  (11) 


I 


•  I 

Recall:  a  ■  1110  ft/sec 


H  «  .414 


! 


I 


Equation  (Q-2)  Is  used  to  convert  the  contours  in  Fig.  11  to 
absolute  vorticity  for  subsequent  calculation  of  secondary  circulation. 
Tabulated  results  are  contained  in  Table  D-I. 


Table  D-I  Absolute  Vorticlty  for  Three  Circumferential  Positions 

on  ufl'  Map 


Radial  I 

Radial 

II 

r'  /  r  (in) 

P' 

ap'/ar'  (sec-1) 

P 

» 

aP'/ar* 

we 

(sec"1) 

.3 

2.43 

.91 

.35  -2554.9 

.9 

.05 

-520 

.325 

2.63 

.93 

-1770.8 

.05 

.35 

2.83 

.94 

-1567.6 

.05 

.375 

3.03 

.95 

-1416.3 

.05 

.4 

3.23 

.96 

-1297.8 

.05 

.425 

3.43 

.97 

-1201.6 

.06 

-624.1 

.45 

3.63 

.98 

-1121.4 

.9 

.075 

-780.1 

.475 

3.83 

.99 

-1053.2 

.91 

.09 

-656.9 

.5 

4.03 

.99 

-1053.2 

.91 

.1 

-729.9 

.525 

4.23 

1.0 

-994.2 

.91 

.11 

-802.9 

.55 

4.43 

1.0 

-994.2 

.92 

.13 

-767.3 

.575 

4.63 

1.01 

-942.5 

.92 

.14 

-826.3 

.6 

4.83 

1.01 

-942.5 

.93 

.15 

-758.9 

.625 

5.03 

1.02 

-896.7 

.93 

.17 

-860.1 

.65 

5.23 

1.03 

•  ' 

-855.7 

.94 

.20 

-895.7 

.675 

5.43 

1.04 

.35  -818.8 

.94 

.23 

-1030.1 

.7 

5.63 

1.05 

.75  -1682.7 

.95 

.25 

-1011.6 

.725 

5.83 

1.06 

1.0  -2156.2 

.97 

.5 

-1716.6 

.75 

6.04 

1.08 

l.i 

2  -2402.2 

1.0 

1.0 

-2840.7 

.775 

6.24 

1.10 

i.; 

2  -2243.5 

1.04 

1.5 

-3509.2 

.8 

6.44 

1.17 

1.0  -1522.4 

1.10 

2.0 

-3739.2 

.825 

6.64 

1.18 

.75  -1112.4 

1.17 

1.5 

-2283.6 

.85 

6.84 

1.19 

0 

0 

1.18 

.5 

-741.6 

.875 

7.04 

1.20 

-.25  352.6 

1.20 

0 

0 

.9 

7.24 

1.18 

-.75  1112.4 

1.18 

-.5 

+741.6 

.925 

7.44 

1.17 

-1.0  1522.4 

1.17 

-1.0 

1522.4 

.95 

7.65 

1.14 

-1.! 

5  2480.2 

1.1 

-3.0 

5608.9 

.975 

7.85 

1.1 

-3.0  5608.9  ! 

1.05 

-3.2 

7179.8 

1.0 

8,05 

Table  D-I  (continued) 


r'  /  r  (in) 

P' 

.3 

2.43 

c 

•  w 

1 

.325 

2.63 

.35 

2.83 

.375 

3.03 

.4 

3.23 

.425 

3.43 

.45 

3.63 

.475 

3.83 

.5 

4.03 

.525 

4.23 

.55 

4.43 

.575 

4.63 

.6 

4.83 

.625 

5.03 

.65 

5.23 

.90 

.675 

5.43 

.90 

.7 

5.63 

.91 

.725 

5.83 

.94 

.75 

6.04 

.97 

.775 

6.24 

1.0 

.8 

6.44 

1.04 

.825 

6.64 

1.17 

.85 

6.84 

1.18 

.875 

7.04 

1.2 

.9 

7.24 

1.18 

.925 

7.44 

1.17 

.95 

7.65 

1.04 

.975 

7.85 

1.0 

1.0 

8.05 

Radial  III 
3P'/3r'  <*>e  (sec 


-.2  +2080.3 


V 

-.2 

» 

+2080.3 

0 

0 

+.3 

-3120.5 

.4 

-4160.7 

.5 

-5200.8 

.5 

-5200.8 

.5 

.5 

.5 

-3649.9 

1.0 

-4478.9 

1.5 

-5149.8 

2.0 

-5681 .4 

2.3 

-5380.8 

1.7 

-2588.0 

0 

0 

-.3 

+423.2 

-.75 

+1112.4 

-1.5 

+2283.6 

-3.0 

+7018.4 

-3.2 

+9090.3 
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/I  sin  23?  -  sin  23i  . 

s  c  U  - Z - 1  -  2c  cos 

m\c  cos  e-, 


«=  55° 


c  *  .3316 


.0480  c„i  ft/sec 


(H-5) 


—  =  S  cos  61  cnl  (tan  -  tan  3j) 


*  .0336 


(E-6) 


Tabulated  results  for  Eqns.  (E-l)  through  E-6)  using  previously 
detc  mined  data  are  contained  in  Tables  E-I,  E-II,  and  E-III. 
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“e  , 
(sec"1) 

C$1, 

(sec"1) 

(sec-1 ) 

5S2 

(sec"1) 

diydr 

(ft/sec) 

dr2/dr 

(ft/sec) 

d^/dr 

(ft/sec) 

-2554.9 

+2364.1 

+1076.3 

+1650 

+121.7 

51.2 

+36.3 

-1770.8 

+1721.8 

+627.0 

+1305.7 

+96.3 

29.8 

+21.1 

-1567.6 

+1555.4 

+510.6 

+1216.5 

+89.7 

24.5 

+17.2 

-1416.3 

+1431.4 

+423.9 

+1150.1 

+84.8 

20.1 

+14.3 

-1297.8 

+1334.4 

+356.0 

+1098.1 

+81.0 

16.9 

+12.0 

-1201.6 

+1255.6 

+300.9 

+1055.9 

+77.9 

14.3 

+10.1 

-1121.4 

+1189.9 

+254.9 

+1020.6 

+75.3 

12.1 

+8.5 

-1053.2 

+1134.0 

+215.8 

+990.7 

+73.1 

10.3 

+7.3 

-1053.2 

+1134.0 

+215.8 

+990.7 

+73.1 

10.3 

+7.3 

-994.2 

+1085.7 

+182.0 

+964.8 

+71.2 

8.6 

+6.14 

-994.2 

+1085.7 

+182.0 

+964.8 

+71.2 

8.6 

+6.14 

-942.5 

+1043.4 

+152.4 

+942.1 

+69.5 

7.2 

+5.12 

-942.0 

+1043.4 

+152.4 

+942.1 

+69.5 

7.2 

+5.12 

-896.7 

+1005.8 

+126.2 

+922.0 

+68.0 

6.0 

+4.2 

-855.7 

+972.3 

+102.7 

+904.0 

+66.7 

4.8 

+3.5 

-818.8 

+942.0 

+81.5 

+887.8 

+65.5 

3.9 

+2.75 

-1682.7 

+1649.6 

.+576.5 

+1267.1 

+93.5 

27.4 

+19.4 

-2156.2 

+2037.5 

+847.9 

+1474.9 

+108.8 

40.3 

+28.7 

-2402.2 

+2239.0 

+988.8 

+1582.9 

+116.8 

47.0 

+33.3 

-2243.5 

+2109.0 

+897.9 

+1513.2 

+111.6 

42.6 

+30.3 

-1522.4 

+1518.3 

+484.7 

+1196.7 

+88.3 

23.0 

+16.3 

-1112.4 

+1182.5 

+249.8 

+1016.7 

+75.0 

11.8 

+8.4 

0 

+271.4 

-387.6 

+528.4 

+38.9 

-18.4 

-13.0 

352.6 

-17.4 

-589.6 

+373.6 

+27.5 

-28.0 

-19.8 

1112.4 

-639.7 

-1025.0 

+40.0 

+2.9 

-48.8 

-34.5 

1522.4 

-975.5 

-1259.9 

-139.9 

-10.3 

-60.0 

-42.5 

2480.2 

-1760.1 

-1808.7 

-560.4 

-41.3 

-86.0 

-61.0 

5608.9 

-4322.8 

-3601.4 

-1933.9 

-142.7 

-171.5 

-121.0 
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nu ~JJl-  u  i.'pi^^'.  v  - 

I  Table  E- 


.85 

.9 

.95 

1.0 


Relative  Vorticity  and  Circulation  Components  at  Radial  II 


(sec-1) 


CS1 

(sec"1 ) 


cnl 

(sec“l) 


5S2 

(sec'1 ) 


dr^/dr  d^/dr 
(ft/sec)  (ft/sec) 


dr3/dr 

(ft/sec) 


-520 

-520 

-520 

-520 

-520 

-624.1 

-780.1 

-656.9 

-729.2 

-802.9 

-767.3 

-826.3 

-758.9 

-860.1 

-895.7 

-1030.1 

-1011.6 

-1716.1 

-2840.7 

-3509.2 

-3738.6 

-2283.6 

-741.6 

0 


+697.3 


+782.3 

+910.3 

+809.4 

+869.2 

+929.0 

+899.8 

+948.2 

+893.0 

+975.9 

+1005.0 

+1115.1 

+1100.0 

+1677.0 

+2598.2 

+3145.7 

+3334.5 

+2141.8 

+878.8 

+271.4 


-89.6 


-29.9 

+59.3 

-11.1 

+30.6 

+72.4 

+52.0 

+85.8 

+47.2 

+105.2 

+125.6 

+202.6 

+192.0 

+595.7 

+1240.1 

+1623.1 

+1755.1 

+920.9 

+37.3 

-387.6 


+756.6 


+802.3 

+870.8 

+816.7 

+848.8 

+880.8 

+865.2 

+891.1 

+861.5 

+905.9 

+921.5 

+980.6 

+972.4 

+1281.7 

+1775.4 

+2068.9 

+2170.0 

+1530.9 

+853.9 

+528.4 


+55.8 


+59.2 

+64.2 

+60.2 

+62.6 

+65.0 

+63.8 

+65.7 

+63.5 

+66.8 

+68.0 

+72.3 

+71.7 

+94.5 

+131.0 

+152.6 

+160.1 

+112.9 

+63.0 

+38.9 


-4.26 

-3.0 

-1.4 

-1.0 

+2.8 

+1.9 

-.5 

-.3 

+1.4 

+1.0 

+3.4 

+2.4 

+2.4 

+1.7 

+4.0 

+2.8 

+2.2 

+1.6 

+5.0 

+3.5 

+5.9 

+4.2 

+9.6 

+6.8 

+9.1 

+6.4 

+28.3 

+20.0 

+59.0 

+41.7 

+77.2 

+54.6 

+83.5 

+59.1 

+43.8 

+31.0 

+1.7 

+1.2 

-18.4 

-13.0 

+741.6 
+1522.4 
+5608.9 
+71 79 i 8 


-336.0 

-975.5 

-4322.8 

-5609.5 


-812.5 

-1259.9 

-3601.4 

-4501.6 


+202.8 

-139.9 

-1933.9 

-2623.5 


+14.9 

-10.3 

-142,7 

-193.6 


-38.6 

-59.9 

-171.4 

-214.2 


-27.3 

-42.4 

-121.3 

-151.7 
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